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ABSTRACT 


The  deflection  behavior  of  composite  beams  under 
service  load  conditions  is  investigated.  The  effects  of 
load-slip  behavior  of  stud  shear  connectors,  shear  defor¬ 
mation,  and  degrees  of  shear  connection  are  considered. 

Basic  equilibrium  equations  in  terms  of  displace¬ 
ment  and  slip  are  formulated  for  composite  beams  with 
partial  shear  connection  and  the  classical  solutions  are 
obtained  for  simply  supported  beams  subjected  to  concentrated 
and  uniformly  distributed  load.  A  numerical  technique, 
based  on  numerical  integration  of  slip  strain,  is  developed 
for  continuous  or  simply  supported  beams  with  a  linear  or 
nonlinear  load-slip  relationship.  This  technique  involves 
a  corrective  iterative  procedure  for  slip  applied  to  an 
assumed  initial  slip,  and  the  numerical  integration  tech¬ 
nique  using  a  fourth  order  Runge-Kutta  method. 

Solutions  obtained  from  the  numerical  analysis 
are  compared  with  theoretical  and  experimental  results  of 
other  investigators.  The  numerical  technique  satisfactorily 
predicts  the  deflection  and  slip  up  to  service  load  con¬ 
ditions.  The  effects  of  partial  and  full  shear  connection, 
linear  and  nonlinear  load-slip  characteristics  of  the  shear 
connectors,  cracking  of  concrete,  and  the  type  of  construction 
method  are  significant  in  the  evaluation  of  deflection  of 
composite  beams. 
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CHAPTER  I 


INTRODUCTION 


1 .  1  Introductory  Remarks 


One  significant  advantage  of  composite 
construction  is  that  beam  deflections  are  considerably 
reduced  as  a  result  of  increased  stiffness.  Composite 
behavior  is  dependent  on  the  shear  connection  between 
the  component  parts. 

1 . 2  Composite  Action 

Shear  connectors  attached  to  the  steel 
section  and  embedded  in  the  concrete  slab  resist  shear 
at  the  interface  of  the  slab  and  steel  section,  thus 
creating  an  interaction  between  two  components.  The  degree 
of  interaction  may  be  complete  or  partial.  A  typical 
composite  beam  is  illustrated  in  Fig.  1.1.  If  there  is  no 
shear  connection  or  friction  at  the  interface  between  the 
concrete  and  steel,  the  compatibility  requirement  reduces 
to  equal  vertical  deflection  of  two  components.  Therefore 
the  applied  load  is  divided  between  two  components  in 
proportion  to  their  stiffnesses  and  each  component  acts  as 
if  it  were  an  isolated  member.  Assuming  the  materials  can 
resist  tensile  stress,  equal  and  opposite  strains  develop 
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at  the  top  and  bottom  fibers  of  each  component ,  as 
illustrated  in  Fig.  1.2(a).  However,  tensile  cracks  may 
occur  in  the  concrete  slab  as  concrete  is  relatively 
weak  in  tension.  The  difference  in  strain  between  the 
adjacent  fibers  of  concrete  and  steel  at  the  interface 
is  known  as  slip  strain,  since  it  results  from  a  horizontal 
slip  that  occurs  between  the  two  components. 

1.2.1  Full  Interaction 

If  the  concrete  slab  and  steel  section  are 
joined  together  by  infinitely  stiff  shear  connectors, 
the  two  members  behave  as  a  unit  with  no  slip  occurring 
at  the  interface.  The  vertical  deflections  of  the 
components  are  equal  at  any  point  along  the  length  of  the 
beam.  Slip  and  slip  strains  are  zero  everywhere  along 
the  interface  and  it  can  be  assumed  that  plane  sections 
remain  plane  during  bending.  This  condition  is  known  as 
full  interaction,  and  is  illustrated  in  Fig.  1.2(c). 

Bending  stresses  and  deflections  due  to  service  loads  for 
the  condition  of  complete  interaction  are  usually  determined 
by  employing  simple  elastic  beam  theory  using  section 
properties  associated  with  the  transformed  cross  section. 

In  practice  the  majority  of  design  methods 
for  composite  beams  are  based  on  the  assumption  of  full 
interaction.  However,  CSA  S16. 1-1974,  "Steel  Structures 
for  Buildings  -  Limit  State  Design",  (  3  )  also  specifies 
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methods  for  composite  design  based  on  partial  interaction. 

1.2.2  Partial  Interaction 

The  assumption  of  full  interaction  is  only 
valid  if  there  is  no  relative  movement  or  slip  at  the 
steel-concrete  interface.  It  is  generally  assumed  that 
the  horizontal  shear  force  at  the  interface  is  transmitted 
only  by  shear  connectors  and  not  by  bond.  Due  to  the 
compressibility  of  the  concrete  slab  and  the  flexibility  of 
the  connectors,  the  shear  force  cannot  be  transmitted 
without  some  slip  occurring,  and  therefore  the  interaction 
must  be  partial,  or  incomplete.  The  strain  distribution 
relating  to  this  type  of  behavior  is  shown  in  Fig.  1.2(b). 
In  every  composite  structure,  interaction  is  less  than 
complete,  irrespective  of  the  relative  strength  of  the 
component  parts.  It  is  therefore  important  to  understand 
how  the  behavior  of  a  composite  beam  is  modified  by  the 
presence  of  slip. 

1 . 3  Deformation  of  Composite  Beams 

Beam  deformation  is  caused  by  bending  and 
shear.  Composite  beam  deformation  is  also  influenced  by 
the  slip  along  the  interface  of  the  two  components. 
Stiffness  rather  than  strength  is  the  governing  criteria 
for  deflection  at  service  loads.  Stiffness  is  influenced 
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by  the  geometry  of  the  cross  section  and  the  material 
properties  of  the  components  involved. 

Bending  deflections  are  dependent  on  the 
moment-curvature  relationship.  Shear  deformation  can  be 
evaluated  independently  of  bending  deformation  by  applying 
the  principle  of  virtual  work.  Deflection  due  to  slip 
depends  on  the  stiffness  of  the  shear  connection  which  is 
defined  by  load-slip  relationship  for  the  type  of  shear 
connector . 

In  continuous  beams,  the  flexural  stiffness 
in  the  negative  moment  region  is  different  from  that  in 
the  positive  moment  region.  In  the  negative  moment  region, 
due  to  tensile  cracking  of  concrete,  the  moment  of  inertia 
of  the  cross-section  is  reduced  considerably  resulting  in  a 
reduction  in  stiffness  which  influences  the  deflection. 
Improper  distribution  of  shear  connectors  or  faulty  con¬ 
nection  may  increase  slip  at  the  interface.  Premature 
yielding  of  shear  connectors  produces  a  non-linear  load- 
slip  relationship  at  working  loads  which  influences  the 
magnitude  of  slip  and  hence  the  deflection. 

1 . 4  Scope 


The  present  study  investigates  the  deflection 
behavior  of  composite  beams  under  service  load  conditions. 

The  effect  of  the  load-slip  behavior  of  stud  shear  connectors 
on  deflection  is  considered.  The  effects  of  partial  and  full 
shear  connection  and  shear  deformation  are  examined.  A 
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numerical  analysis  technique  is  developed  for  the  analysis 
of  single  span  and  continuous  composite  beams.  A  further 
objective  is  to  compare  deflections  produced  in  shored 
and  unshored  construction. 
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CHAPTER  II 


REVIEW  OF  PREVIOUS  RESEARCH 

2 . 1  Load-Slip  Relationship  for  Shear  Connectors 

2.1.1  Push-Out  Test 

The  push-out  test  has  been  used  at  least  since 
1930  (26)  as  a  means  of  evaluating  the  load  carrying 
capacity  of  shear  connectors  and  the  load-slip  relationship 
for  shear  connectors  is  normally  determined  from  such  a 
test.  Fig.  2.1  shows  a  typical  push-out  test  specimen 
and  Fig.  2.2  illustrates  typical  results  obtained  from 
such  a  test . 

In  push-out  tests,  one  approach  considers 
that  the  function  of  shear  connectors  is  to  control  the 
magnitude  of  slip  between  the  concrete  slab  and  the  steel 
section.  The  load  carried  by  the  connectors  at  some 
limiting  slip  is  defined  as  the  "useful  load  capacity" 

(25).  A  second  approach  considers  that  the  function  of 
shear  connectors  between  sections  of  zero  and  maximum 
moment  is  to  transfer  across  the  concrete-steel  interface 
the  maximum  compressive  force  that  can  be  developed  in 
the  concrete  slab,  or  the  maximum  tension  force  in  the 
steel  section  without  consideration  of  the  magnitude  of 
slip  which  has  occurred. 
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2.1.2  Linear  Load-Slip  Behavior 


Newmark  et  al  (17),  in  interpreting  their 
push-out  test  results,  concluded  that  the  amount  of  slip 


modulus,  K,  as 


K  = 
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Q 


Y 

q£ 
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1  app 
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2.1.3  Non-Linear  Load-Slip  Behavior 


On  the 

tests,  Slutter  and 
load  capacity  of  a 


basis  of  a  study  of  numerous  push-out 
Driscoll  (22)  defined  the  ultimate 
stud  connector  as 


where  Q  is 
u 

connect  or , 

and  U  t  is 
cyl 


0  =  930d2  /U  “ 

'u  cyl 

the  ultimate  load  capacity  of  a  shear 
d  is  the  diameter  of  the  connector  in 
the  28-day  concrete  cylinder  strength 


(2.3) 
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Van  Dalen  (24)  proposed  the  following 
relationship  between  slip,  applied  shear  force  and 
foundation  modulus: 


n-2  nD-0-75 

Y  =  3.7  x  10  QR 


(2.4) 


where  Q  is  the  applied  shear  force  in  pounds,  and  R  is 
foundation  modulus  in  psi.  The  foundation  modulus  was 
defined  as  the  modulus  of  elasticity  of  concrete.  Fig 
2.3  illustrates  the  relationship  proposed  by  Van  Dalen 
From  his  tests  he  concluded  that : 


1.  The  ultimate  load  carrying  capacity  of  a 
stud  connector  embedded  in  a  concrete  slab 
is  primarily  dependent  upon  the  tension  in 
the  connector  resulting  from  a  separation  of 
concrete  and  steel  at  the  interface. 

2.  The  ultimate  load  of  a  connector  in  a 
push-out  test  should  not  be  adopted  as  the 
ultimate  load  for  a  connector  in  a  beam  unless 
the  expected  separation  at  the  concrete-steel 
interface  in  the  beam  is  less  than  that  at 
failure  in  a  push-out  test. 


Yam  and  Chapman  (30)  proposed  an  exponential 
relationship  between  load  and  slip  for  a  stud  shear 
connector : 

Q  =  a(l  -  e~bY) 


(2.5) 


where  a  and  b  are  constants.  The  exponential  relationship 
was  obtained  by  establishing  the  best  fit  curve  for 
experimental  results.  By  selecting  two  slip  values  in 
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an  experiment 
is  two  times 
evaluated  as 


al  load-slip  plot  such  that  one  slip 
the  other,  the  constants  a  and  b  can 


Q- 


2Qp  —  Q2 


b  =  —  &n 


Q 


1 


Y  1  Q2  -  Q1 


value 

be 


(2.6) 

(2.7) 


Fig.  2.4  illustrates  this  exponential  load-slip  relationship 
and  defines  values  of  Y^,  Y^  ,  and  •  The  ultimate 
capacity  of  a  shear  connector  was  defined  by  Yam  and 
Chapman  as  the  load  at  which  the  slip  attains  a  limiting 
value  of  0.055  inches. 


2 . 2  Section  Properties 
2.2.1  Positive  Moment  Region 

Properties  of  composite  beams  are  different 
in  positive  and  negative  moment  regions.  In  a  positive 
moment  region  as  shown  in  Fig.  2.5(a),  the  concrete  slab, 
or  a  portion  thereof,  depending  upon  the  location  of 
neutral  axis,  is  in  compression  and  contributes 
significantly  to  moment  resistance.  In  evaluating  the 
elastic  stiffness  of  a  composite  section,  it  is  customary 
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to  transform  the  concrete  slab  area  into  an  equivalent 
area  of  steel  by  applying  the  modular  ratio  (Ec/Es),  as 
a  multiplier  of  the  effective  slab  width. 

Slutter  and  Driscoll  (22)  proposed  that  the 
load  capacity  of  a  connector  in  a  positive  moment  region 
should  be  the  ultimate  load  as  obtained  in  a  push-out 
test ,  whereas  Chapman  and  Balkrishnan  (  5  )  propose  80% 
of  ultimate  capacity. 

2.2.2  Negative  Moment  Region 

In  the  vicinity  of  an  interior  support,  a 
continuous  composite  beam  is  subjected  to  negative  bending 
moment  which  produces  tension  in  the  slab.  If  the  concrete 
cracks,  the  composite  section  consists  of  the  longitudinal 
slab  reinforcement  and  the  steel  section  as  shown  in 
Fig.  2.5(b).  The  cracked  concrete  acts  as  a  medium  for 
transferring  horizontal  shear  forces  required  to  develop 
tension  stress  in  the  longitudinal  steel.  Thus  the 
stiffness  is  significantly  lower  than  that  in  a  positive 
moment  region.  This  decreased  stiffness  may  result  in 
an  increase  of  as  much  as  25%  in  the  value  of  the  positive 
moment  based  on  uniform  stiffness.  The  elastic  stress  and 
strain  distribution  for  composite  beams  in  positive  and 
negative  moment  region  is  illustrated  in  Fig.  2.6. 
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Siess  (21)  reported  results  of  tests  on 
two-span  continuous  bridge  beams,  one  with  shear 
connectors  provided  along  the  total  length  of  the  beam, 
and  the  other  with  shear  connectors  only  in  the  region 
of  positive  moment.  Longitudinal  slab  reinforcement 
was  provided  in  the  negative  moment  region  of  both  beams. 
Siess  concluded  that,  at  ultimate  moment,  complete 
interaction  existed  in  the  positive  moment  region  of  the 
beam  with  shear  connectors  throughout  the  length,  whereas 
partial  interaction  was  present  in  the  other  beam.  The 
tests  thus  showed  that  shear  connectors  are  effective  in 
the  negative  moment  region  and  they  are  necessary  to 
achieve  effective  composite  action. 

Van  Dalen  (24)  concluded  that  stud  shear 
connectors  form  a  satisfactory  shear  connection  in  a 
negative  bending  moment  region.  However,  he  concluded, 
their  capacity  is  somewhat  less  than  that  attained  in 
a  positive  moment  region. 

Wastlund  and  Ostlund  (29)  tested  composite 
beams  loaded  so  that  the  concrete  slab  was  on  the  tension 
side.  One  beam  had  only  two  bow-shaped  shear  connectors, 
while  a  second  beam  contained  no  shear  connectors  but 
contained  slab  reinforcement  bent  down  and  welded  to  the 
steel  beam  at  the  ends.  The  third  beam  had  a  prestressed 
concrete  slab  with  three  bow-shaped  connectors  at  each 


end.  The  authors  concluded  that: 
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1.  In  a  composite  beam  subjected  to  negative 
bending  moment ,  the  concrete  cannot  be 
assumed  to  act  compositely  with  the  steel 
section,  as  the  concrete  cracks  at  very 
small  tensile  stresses. 

2.  If  the  connection  between  the  steel  section 
and  slab  is  sufficient  (there  is  no  indication 
of  what  might  be  a  sufficient  connection)  the 
longitudinal  reinforcement  acts  jointly  with 
the  steel  section. 

3.  The  ultimate  load  capacity  of  shear  connectors 
in  a  negative  moment  region  is  considerably 
smaller  than  that  found  for  similar  shear 
connectors  in  positive  moment  regions. 

The  amount  of  longitudinal  reinforcement 
affects  the  behavior  of  composite  beams  in  the  negative 
moment  region.  Piepgrass  (20)  observed  that  the  ratio  of 
experimental  to  theoretical  ultimate  moment  decreases  with 
an  increase  of  amount  of  longitudinal  reinforcement  in 
negative  moment  regions.  Davison  (  6  )  concluded  that  the 
amount  of  longitudinal  reinforcement  increases  the 
negative  moment  capacity  but  reduces  the  rotation  capacity 
of  plastic  hinges  in  the  negative  moment  region.  In  his 
study  the  amount  of  longitudinal  reinforcement  was  0,111 


to  0.232  times  the  area  of  the  steel  section. 
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2 . 3  Analytical  Studies  of  Simply  Supported  Composite 


Beams 


2.3.1  Elastic  Analysis 


Newmark  et  al  (17)  developed  a  closed  form 


solution  for  a  simply  supported  composite  beam.  They 
assumed  a  linear  stress-strain  relationship  for  steel  and 
concrete  and  a  linear  load-slip  relationship  for  shear 
connectors.  Their  analysis  is  based  on  the  following 
assumptions : 


1.  The  shear  connection  between  slab  and  the 
beam  is  assumed  to  be  continuous  along  the 
length  of  the  beam. 

2.  The  amount  of  slip  occurring  in  the  shear 
connector  is  assumed  to  be  directly  propor¬ 
tional  to  the  load  transmitted. 

3.  The  distribution  of  strain  throughout  the 
depth  of  slab  and  steel  section  is  assumed 
linear . 

4.  The  steel  section  and  the  concrete  slab  are 
assumed  to  deflect  equal  amounts  at  all 
points  along  the  span  length. 


They  defined  slip,  or  the  relative  movement  between  the 
concrete  slab  and  steel  section  interface,  as 


q  £  _  £ 


dF 


(2.8) 


K  K 


dx 


•» 
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where  F  is  the  compressive  force  in 
The  rate  of  change  of  slip  is  equal 
between  strain  in  the  concrete  slab 
at  the  interface. 


the  concrete  slab, 
to  the  difference 
and  the  steel  section 


T.  =  £  _  £ 

dx  b  s 


(2.9) 


where  £  is  the  strain  in  the  steel  and  £  is  the  strain 
b  s 

in  the  concrete. 

Using  a  linear  stress  relationship  they 
developed  the  following  second  order  differential  equation 
for  force  in  terms  of  moment ,  shear  connector  modulus  and 
section  properties: 

F"  +  -j-  (qM  -  f2F)  =  0  (2.10) 

where  M  is  bending  moment ,  and 
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s 

T 
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moment 

of 

slab  and  beam,  respect 

ively , 

and 

y 

is  the 

between  centroid  of  the  concrete  slab  and  the 
the  steel  section. 

Newmark  et_  al_  (17 )  also  presented 
for  Equation  2.10  for  a  simply  supported  beam 


(2. lib) 

y  of  concrete 
of  intertia 
distance 
centroid  of 

a  solution 
subjected  to 


a  single  concentrated  load  only.  They  found  close 


agreement  between  theoretical  results  and  experimental 
results  obtained  from  a  series  of  tests. 
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2.3.2  Inelastic  Analysis 


Yam  and  Chapman  (30)  studied  ultimate  load 
behavior  of  simply  supported  composite  beam  using  a 
non-linear  stress-strain  relationship  and  inelastic 
load-slip  behavior  for  the  shear  connectors.  Based  on 
the  equilibrium  of  a  segment  of  the  slab,  the  shear  flow, 
q,  at  the  interface  of  the  concrete  slab  and  steel  section 
may  be  defined  as 


where  C  is 
The  strain 
horizontal 


Q  = 


dC 

dx 


the  compressive  force  as  shown  in 

difference  at  the  interface  due  t 

movement  of  two  components  may  be 

P  =  dy 
d  dx 


(2.12) 

Fig.  2.7. 
o  the  relative 
expressed  as 

(2.13) 


Using  an  exponential  load-slip  relationship  as  described 
in  Section  2.13,  they  derived  the  following  two  first 
order  simultaneous  equations  with  the  dependent  variables 
C  and  Y  as 


C' 


dC 

dx 


e 


-by 


) 


(2.14) 


Y  ' 


dY 

dx" 


f  ( M ,  C ) 


(2.15) 


where  f(M,C)  is  an  implicity  function  of  moment  and 
compressive  force. 
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By  eliminating  y,  the  following  second  order 
differential  equation,  with  C  as  a  dependent  variable, 
results : 


C"  +  C' 


d£ 


i  dx 


-  b  f  ( M ,  C ) 


+ 


ab 


i 


f ( M , C )  =  0  (2.16) 


This  equation  may  be  solved  by  the  predictor-corrector 
method  of  numerical  integration.  For  a  simply  supported 
beam  a  value  of  slip  is  assumed  at  one  end  and  force  is 
computed  at  the  other  end.  If  the  computed  force  is  zero, 
then  the  assumed  slip  value  is  correct,  otherwise 
successive  corrections  are  applied.  Yam  and  Chapman  found 
good  correlation  between  theoretical  and  experimental 
results  up  to  98%  of  the  failure  load,  but  there  was 
significant  difference  at  connector  failure. 


2 . 4  Analytical  Studies  of  Continuous  Composite  Beams 
2.4.1  Elastic  Analysis 

Plum  and  Horne  (19)  developed  a  closed  form 
solution  for  a  continuous  composite  beam  of  two  unequal 
spans  with  equal  concentrated  loads  on  each  span.  They 
used  a  linear  load-slip  relationship  for  the  shear 
connectors,  elastic  properties  for  the  concrete  slab  and 
steel  section,  and  partial  interaction  between  the  concrete 
slab  and  steel  section.  They  derived  the  following 
governing  fourth  order  differential  equation,  in  terms 
of  deflection  as  the  dependent  variable: 
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v 


IV 


kv" 


w 


yM 


(2.17) 


(SEI) 


EA(EEI ) 


(2.18) 


EA(ZEI ) 


where  y  is  the  interface  stiffness  (shear  connector 
modulus),  v  is  the  deflection,  M  is  the  bending  moment, 
w  is  the  uniformly  distributed  load,  (EEI)  and  EA  are 
composite  section  propeties,  and  z  is  the  distance  between 
the  centroids  of  the  concrete  slab  and  steel  section. 


Beams  with  constant  and  variable  flexural 


stiffness  were  studied.  For  constant  flexural  stiffness, 
the  cracking  of  the  concrete  slab  in  the  negative  bending 
moment  region  was  ignored  and  the  composite  beam  was  assumed 
to  be  of  uniform  cross  section  over  the  entire  length. 
Variable  stiffness  was  treated  by  means  of  an  equivalent 
haunched  beam  in  which  depth  of  the  beam  along  the  length 
varied  in  accordance  with  the  moment  of  inertia  of  the 
actual  beam  section.  In  other  words,  the  variable  stiffness 
did  not  arise  from  cracking  of  concrete  in  the  negative 
moment  region.  A  comparison  of  theoretical  and  experimental 
deflection  values  indicated  good  agreement. 

2.4.2  Inelastic  Analysis 


Yam  and  Chapman  (31)  extended  their  analysis 


of  a  simply  supported  beam  to  a  two-span  symmetric 
continuous  beam  which  could  be  modeled  as  a  single  span 
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propped  cantilever.  Equal  concentrated  loads  were  applied 
to  each  span  and  the  same  composite  section  was  used  as 
for  the  simple  span  case.  The  effect  of  shear  connector 
spacing  and  types  of  loading  on  deflection  and  slip  were 
studied.  They  used  the  governing  equation  2.16  and  the 
same  integration  procedure  as  for  the  simply  supported  beam. 
In  the  solution  for  a  continuous  beam,  values  of  both  slip 
at  the  exterior  support  and  bending  moment  at  the  interior 
support  must  be  initially  assumed.  Therefore,  more 
iterations  requiring  considerable  computational  time  must 
be  performed  to  obtain  correct  values. 

Hamada  and  Longworth  (9)  suggested  the  com¬ 
putational  time  required  in  Yam  and  Chapman's  procedure 
can  be  significantly  reduced  by  assuming  that  the  slip 
strain  is  constant  along  the  shear  span.  This  assumption 
is  based  on  linear  slip  distribution  and  is  satisfactory 
for  a  simply  supported  beam  but  is  not  valid  for  a 
continuous  beam. 

2 . 5  Effect  of  Slip 

The  slip  at  the  interface  of  the  slab  and  the 
steel  section  affects  stress,  strain  and  deflection  at  all 
sections  along  the  beam.  According  to  Siess  (21)  the 
effect  of  slip  on  strain  is  a  maximum  at  the  interface  and 
minimum  at  the  bottom  of  the  steel  section  where  the  strain 
is  maximum.  His  tests  indicated  that  the  effect  of  slip 
on  the  distribution  of  strain  was  a  relatively  localized 
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effect  confined  to  the  region  of  the  applied  concentrated 
load . 

Johnson  (11)  reported  that  within  the  elastic 
range,  slip  may  change  the  stress  distribution  by  as  much 
as  5%  and  may  increase  deflection  by  as  much  as  13%. 

Newmark  et_  al  (17)  obtained  similar  results  in  full  scale 
beam  tests. 

Plum  and  Horne  (19)  stated  that  incomplete 
interaction  may  increase  the  deflection  as  much  as  40%. 

Their  test  results  indicate  that  slip  may  increase 
deflection  as  much  as  50%  and  decrease  the  compressive  force 
in  the  slab  up  to  20%. 

Hamada  and  Longworth  (  9  )  proposed  an  analysis 
for  computing  deflections  of  continuous  composite  beams 
which  included  the  effect  of  shear  and  slip.  They  found 
that  shear  deformations  were  significant  and  their 
analytical  results  were  in  close  agreement  with  actual 


deflections  determined  in  tests. 


* 
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FIGURE  2.1  PUSH-OUT  TEST  SPECIMEN 
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FIGURE  2.2  LOAD-SLIP  RELATIONSHIP  FROM  PUSH-OUT  TEST 

(CHAPMAN  AND  BALAKRISHNAN ) 
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FIGURE  2.3  LOAD-SLIP  RELATIONSHIP  (VAN  DALEN) 
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FIGURE  2.4  LOAD-SLIP  RELATIONSHIP  (YAM  AND  CHAPMAN) 
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FIGURE  2.6  ELASTIC  STRAIN  AND  STRESS  DISTRIBUTION 


FOR  COMPOSITE  BEAM 
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(a)  STRESS  &  STRAIN  ACROSS  A  SECTION 
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ELEMENT  OF  SLAB 


(c)  ELEMENT  OF  COMPOSITE  BEAM 


FIGURE  2.7 


CONDITIONS  IN  AN  ELEMENT  OF  A  COMPOSITE  BEAM 


CHAPTER  III 


FORMULATION  AND  CLASSICAL  SOLUTION 
OF  BASIC  EQUATIONS  FOR  COMPOSITE  BEAMS 

3 . 1  Formulation  of  Basic  Equations 

3.1.1  Assumptions  for  Displacements  and  Stresses 


In  the  analysis,  the  following  assumptions 

are  made : 

1.  The  distribution  of  strain  is  linear  over  the 
depth  of  slab  and  the  depth  of  steel  beam, 
respectively . 

2.  The  shear  connection  between  the  slab  and 
steel  beam  acts  as  a  continuous  medium  along 
the  length  of  the  beam. 

3.  Concrete  has  no  tensile  strength. 

4.  It  is  assumed  that  the  reinforcement  bars  are 
placed  at  one  depth  in  the  concrete  slab. 

5.  The  stress/strain  curves  for  steel  are  the  same 
in  tension  and  in  compression. 

6.  Within  service  load  range  the  stress/strain  for 
concrete  and  steel  are  linear. 

7.  The  concrete  slab  and  steel  beam  deflect  equally 
at  all  points  along  the  beam  so  that  at  any 
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cross  section  they  have  equal  curvature 
and  that  uplift  forces  are  resisted  by  the 
shear  connectors  without  separation  and  do 
not  affect  the  behavior  of  composite  beams. 

Fig.  3.1  shows  a  portion  of  a  composite  system 
with  the  slab  spanning  over  several  equally  spaced  beams. 

In  the  transverse  direction  the  slab  is  considered  to  act 
as  a  continuous  one-way  slab  supported  by  the  beams.  A 
portion  of  the  slab,  of  width  be,  acts  compositely  with 
each  steel  section.  The  coordinate  system  and  reference 
dimensions  for  a  typical  composite  beam  are  shown  in 
Fig.  3.2(a)  and  the  assumed  deformation  and  displacements 
are  shown  in  Fig.  3.2(b). 

As  a  result  of  partial  interaction,  cross 
section  ABCD  which  is  plane  before  deformation  will  not 
remain  plane  after  deformation.  The  slip,  s,  creates  a 
discontinity  at  the  interface  of  the  concrete  and  steel 
as  shown  in  Fig.  3.2(b). 

The  strains  at  any  point  in  the  beam  may  be 
determined  directly  from  the  horizontal  displacement  u 
and  the  vertical  displacement  v  through  the  use  of 
standard  beam  assumptions  and  strain  displacement  equations. 
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Assuming  plane  sections  remain  plane,  the  horizontal 
displacements  in  the  steel  may  be  expressed  as 

u  =  uQ  +  v’y  0<y<  d  (3.1) 

and  the  horizontal  displacements  in  the  concrete  may  be 
expressed  as 

u  =  Uq  +  v'y  +  s  d<y<(d  +  t)  (3.2) 

Combining  Equations  3.1  and  3.2,  and  using  a  step  function 
(  7 )  so  that  the  last  term  of  Equation  3.2  is  zero  in 
the  range  0<y<d,  the  displacement  across  the  entire 
section  can  be  expressed  as 


u  = 

uo 

+  v'y  + 

<y 

,  0 

-  d  >  s 

0  <y  < 

(d  +  t) 

(3.3) 

where 

<y  - 

,  0  , 
d>  has 

a 

value  of  0 

f  or 

a  negat 

i  ve 

argument 

and  1 

f  or 

a  positi 

ve 

argument . 

The 

strain 

in  t 

he 

horizont  al 

direct i 

on 

can  now  be 

dete 

rmined 

f  rom 

the  strain 

displacement  equation  as 

e  =  u '  (3.4) 

x 

Differentiating  Equation  3.3,  yields 

e  =  u ' q  +  v"y  + <  y  -  d>°s '  (3.5) 

Defining  the  reference  axis  strain, 

u ' 0  =  £0 


(3.5a) 
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the  curvature  as 

v"  =  -  <j)  (3.5b) 

and  slip  strain  as 

Srs-  (3.5c) 

Equation  3.5  may  be  written  as 

ex  =  e0  “  +  <y  ~  d>°  s'  (3.6) 

For  a  linear  elastic  stress-strain  relationship,  stress 
is  expressed  as 


a  =  Ee  (3.7) 

x 

Substituting  Equation  3.6  into  Equation  3.7  gives  stress 
in  the  terms  of  displacement  derivatives  as 

a  =  E  (  £q  -  4>y  +  <y  -  d>°  s'  )  (3.8) 


3.1.2  Equilibrium  Equations  in  Terms  of  Displacements 


Three  basic  equilibrium  equations  may  be 
written  for  the  composite  beam  shown  in  Fig.  3.4.  These 
are 

Jf^adA  =0  (3.9) 


-  M 


d 

dx 


■  d 

/  adA 

-  0 


dP 

dx 


Ks 


q 


(3. 10) 


(3.11) 
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where  q  is  the  shear  flow  on  the  interface  between  steel 
and  concrete.  Equation  3.11  is  valid  only  if  there  is  a 
relationship  between  shear  flow  and  slip.  Substituting 
Equation  3.8  into  Equations  3.9  to  3.11,  and  carrying 
out  the  integration,  we  obtain 


E( A  +  A  )  en  -  E ( A  Y  +  A  Y  )tf>  +  EA  s'  =  0  (3.12) 

v  s  c  0  v  c  c  s  s  'v  c  v  ' 


E ( A  Y  +  A  Y  )e_  -  E(I  +  A  Y  2  +  I  +  A  Y  2 )  6 
c  c  s  s  °0  c  cc  s  ss^ 


+  EA  Y  s’  =  -M  (3.13) 

c  c 

EA  e’  -  EA  Y  a'  =  -  Ks  (3.14) 

s^  0  s  sv 

in  which  Y  and  Y  indicate  the  distances  from  neutral 
s  c 

axis  of  the  composite  section  to  the  centroids  of  the 

area  of  steel,  A  ,  and  the  transformed  area  of  concrete, 

s 

Ac  (neglecting  concrete  area  in  tension),  respectively, 

E  is  the  modulus  of  elasticity  of  the  steel,  and  n 
(=  Ec/Eg)  is  the  inverse  modular  ratio.  Equations  3.12 
to  3.14  may  be  simplified  by  shifting  the  x-axis  to  the 
centroid  of  the  transformed  section,  in  which  case  the 
following  definitions  apply: 


(A  Y 
s  s 


A  Y  )  =  0 
c  c ' 


(3.14a) 


and 


(I. 


+ 


A  Y  ‘ 
c  c 


+  I 


A  Y  ‘ 
s  s 


)  =  I 


t 


(3.14b) 


where 

steel 


I  is  the  moment  of  inertia  of  the  steel  about  the 
s 

section  centroidal  axis,  I  is  the  moment  of 

c 
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inertia  of  the  transformed  area  of  the  concrete  slab 


about  its  centroidal  axis  in  compression  and  I  is  the 
moment  of  inertia  of  the  transformed  composite  section 
about  the  centroidal  axis  of  the  composite  section.  Thus 
e q  becomes  the  strain  at  the  composite  section  centroid 


level,  and 


A  +  A  =  A , 
set 


(3.14c) 


is  the  total  area  of  transformed  section.  The  simplified 
forms  of  Equations  3.12  to  3.14  are 


EA,e„  +  EA  s '  =  0 
tO  c 


-  EIJ  +  EA  Y  s'  =  -  M 
t  y  c  c 


EA  e'n  -  EA  Y  f  =  -Ks 
s  0  s  sY 


(3.15) 

(3.16) 

(3.17) 


Equations  3.15  to  3.17  express  the  basic  equilibrium 

conditions,  i.e.,  Equations  3.9  to  3.11,  in  terms  of  the 

du. 


0 


0 


dx 


,  (p  -  - 


, 2 

d  v 


dx 


2 


,  and 


three  displacement  gradients  z 
ds 

s'  =  ;  the  section  properties  of  the  transformed  area 

and  the  stiffness  of  the  shear  connectors,  K.  However, 
in  order  to  obtain  a  solution,  it  is  desirable  to  eliminate 
all  but  one  of  the  displacements  from  the  equations  and 
derive  governing  equations  to  enable  the  solution  of  one 
displacement  independently  from  the  others. 


V 
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3.1.3  Governing  Equation  for  Slip 

By  differentiating  Equations  3.15  and  3.16, 
e'  o  and  cjd  may  be  expressed  in  terms  of  derivatives  of  s 


Thus , 


A 


0 


0 


A 


s  ' 


t 


A 


A 


A  Y  s’  .. 
(f>  =  -CC _  +  M 


t 


El 


A  Y  s" 

a  ,  _  c  c  ,  V 

4>' - - -  + 


t 


El 


(3.18a) 


(3.18b) 


(3.19a) 


(3.19b) 


Substituting  these  expressions  for  e  and  <j)'  into 
Equation  3.17  and  grouping  the  terms  with  a  common  order 
of  derivative  of  slip,  s,  yields 


s”  - 


K 

V 

A  Y 
s  s 

i 

A 

A  Y  Y 

S  E 

I ,  A 

A 

A  Y  Y 

E  A  S 

J.  s  s  c 

t  c 

s  + 

s  s  c 

h  c  A 

T 

A , 

L 

L  t 

t  J 

L  t 

t  J 

Letting 


K 


A 


A  Y  Y 


EA  ^  +  S  S  C 


c  A 


t 


I 


t 


=  a 


and 


A  Y 
s  s 

El  ,  A 
t  c 


1 


A  A  Y  Y 
s  +  s  s  c 


A 


t 


I 


=  6 


(3.20) 


(3.21) 


(3.22) 


Equation  3.20  becomes 


s"  -  a  s  = 


-  6  V 


(3.23) 
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This  is  a  second  order  differential  equation  and  is  the 
governing  equation  for  slip.  When  the  beam  is  statically 
determinate,  V  is  a  known  function  of  x,  and  s  can  be 
determined  by  solving  the  differential  equation.  With 
the  slip  known,  the  curvature  can  be  obtained  directly 
from  Equation  3.15.  Integrating  Equation  3.16  twice  yields 
the  value  of  transverse  deflection  v,  and  integrating 
Equation  3.15  once  yields  the  value  of  Uq.  Thus  the 
complete  solution  may  be  obtained  from  Equations  3.23, 

3.15,  and  3.16. 


3.1.4  Governing  Equation  for  Deflection 


An  alternative  solution  technique  consists 
of  deriving  a  governing  fourth  order  equation  in  terms  of 
the  transverse  displacement  v.  Differentiating  Equation 
3.17  and  expressing  terms  as  displacement  gradients 
yields 

EA  u'"  +  EA  Y  v1V  +  Ks'  =0  (3.24) 

s  0  s  s 


From  Equation  3.16, 


s  T  = 


-M 


EA  Y 
c  c 


A  Y 
c  c 


V' 


(3.25) 


By  substituting  Equation  3.25  in  Equation  3.18a 
and  differentiating  twice,  Uq'"  may  be  obtained  in  terms 
of  v.  Thus 

A  I  v" 

,  _  c  ,  M  t  v 

U  0  A,  lEA  Y  A  Y  ; 

t  c  c  c  c 


(3.26) 
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u 


0 


-w 

EY  A , 
c  t 


+ 


Itv 


IV 


AY 
t  c 


(3.27) 


Substituting  Equations  3.25  and  3.27  into  Equation  3.24 
results  in 


■  I,  +  A  Y  Y  T  A  A 

-v1VEA  A  (  —  r - — -  C)  +  KI.  v"  =  ~  w 

s  c  A  '  t  E  A 

^  L 


Letting  I,  +  A,Y  Y 

EA  A  (  t  A  1  S  C)  =  * 
s  c  v  At  ' 


(3.28) 

(3.29) 


Equation  3.28  becomes 

KI 


v 


IV 


t 


4^ 


V"  = 


-KM 

E  ip 


A  A 
s  c 

V 


w 


(3.30) 


This  is  a  fourth  order  differential  equation 
in  terms  of  deflection.  If  the  beam  is  statically 
determinate,  M  is  a  known  function  of  x  and  w  (the 
uniform  loading)  and  v  can  be  determined  as  the  solution 
to  the  differential  equation.  With  deflection  known,  the 
curvature  may  be  determined  by  differentiation  and  the 
slip  strain  determined  from  Equation  3.25  by  substitution. 
The  slip  may  be  evaluated  by  integrating  Equation  3.25 
and  Uq  is  determined  from  Equation  3.26.  The  complete 
solution  is  therefore  available  from  Equations  3.30, 


3.25,  and  3.26. 
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3 . 2  Slip  Strain  at  Inflection  Point 

Figure  3.5  illustrates  the  neutral  axis 
positions  in  a  continuous  composite  section  under  positive 
and  negative  bending  for  uniform  and  discontinuous  section 
properties.  In  the  negative  bending  moment  region, 
assuming  the  concrete  has  no  tension  strength,  the  concrete 
slab  is  not  effective.  Theoretically  there  is  an  abrupt 
change  in  the  position  of  neutral  axis  at  the  point  of 
inflection  as  shown  in  Fig.  3.5(b).  The  slip  strain  is 
affected  by  this  abrupt  change  in  the  neutral  axis 
position  and  the  abrupt  change  in  effective  area  of 
concrete  slab. 

In  order  to  maintain  continuity  of  the  force 
F  between  positive  and  negative  moment  regions,  the  strain 
at  the  centroid  of  steel  section  must  be  equal  on  both 
sides  of  the  point  of  inflection. 


esL  esR 


(3.31) 


where  subscripts  L  and  R  refer  to  locations  to  the  left 
and  right  of  inflection  point.  In  the  positive  moment 
region  the  strain  at  the  centroid  of  the  steel  section 
can  be  expressed  in  terms  of  the  curvature  and  the  strain 
at  the  neutral  axis  of  the  transformed  section.  From 


Equation  3.6, 


esL  e0L 


L  sL 


(3.32a) 
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Similarly  for  the  negative  moment  region, 


EsR  £0R  VsR 


(3.32b) 


Now  substituting  values  of  £qL>  £0R’  anc*  ^R  from 
Equations  3.18(a),  and  3.19(a)  into  Equations  3.32(a) 
and  3.32(b),  and  equating  the  results  as  indicated  by 
Equation  3.31,  the  following  relationship  is  obtained: 


A 


Y  tA  tY  t 
cL  sL  cL  cL 


tL 

A 


cR 


tR 


+ 


tL  / 


Y  tA  dY 
sL  cR  cR 

1 1 R 


w 


MY 


L  +  El 


sL 


tL 

MY 


R 


+ 


sR 


El 


(3.33) 


tR 


If  uniform  or  constant  section  properties 
are  assumed  over  the  entire  length  of  beam,  there  will  be 
no  change  in  the  concrete  area  at  the  inflection  point 
and  the  neutral  axis  will  not  shift  its  position.  Hence 
the  section  properties  to  the  left  of  the  inflection 
point  will  be  the  same  as  section  properties  to  the  right. 
From  Equation  3.33  it  is  evident  that  with  uniform  section 
properties  the  slip  strain  to  the  left  of  the  inflection 
point  is  equal  to  that  to  the  right.  Thus  it  can  be 
proved  that  continuity  of  slip  strain  exists  with  constant 
section  properties. 

If  there  is  a  change  in  section  properties  the 
slip  strain  to  the  right  of  the  inflection  point  will  be 
related  to  slip  strain  to  the  left  but  there  will  be 
discontinuity  at  the  inflection  point.  Let 
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A 


lcL  +  YsL^cL^cL 


tL 


tL 


and 


A  D  Y  A  dY  _ 
cR  +  sR  cR  cR 


tR 


tR 


AL 


X 


R 


(3.34) 


(3.35) 


By  substituting  Equations  3.34  and  3.35  into  Equation 
3.33  and  eliminating  moment  terms  (moment  is  zero  at  the 
inflection  point),  the  slip  strain  to  the  right  of  the 
inflection  point  may  be  written  as 


s ' 


R 


X 


s  ' 


R 


L 


(3.36) 


3 . 3  Boundary  Conditions 

In  order  to  evaluate  the  constants  of 
integration  involved  in  the  solution,  boundary  conditions 
must  be  established  for  any  particular  problem.  For 
example,  for  a  simple  span  beam  the  bending  moment, 
stresses  and  strains  are  zero,  i.e. ,  M  =  0,  a  =  0, 

=  0  at  x  =  0  and  at  x  =  L.  Hence,  from  Equation  3.15, 
s'  must  be  zero  and,  from  Equation  3.16,  it  may  be 
concluded  that  <p  is  also  zero  at  x  =  0  and  at  x  =  L. 

From  compatibility  requirements,  it  is 
obvious  that  slip  is  continuous.  Also,  as  explained  in 
Section  3.2,  the  strain  at  the  centroid  of  the  steel 
section  is  continuous  since  the  force  P  in  the  slab  must 
be  continuous  and  therefore  s’  is  continuous  wherever 
section  properties  are  continuous.  For  equilibrium,  the 
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summation  of  the  axial  force  over  the  entire  section 
must  be  zero.  The  boundary  conditions  for  a  simply 
supported  beam  and  a  continuous  beam  are  shown  in 
Fig.  3.6. 


3.4  Linear  Closed  Form  Solutions 


The  governing  equation  for  slip,  Equation 


3.23,  is 


sM  -  a  s  =  -3  V 


(3.37) 


The  particular  solution  of  this  equation  for  a  uniformly 
distributed  load,  w,  is 


P 


3  ,wL  v 

-~2  -  wx) 
a 


(3.38) 


The  homogeneous  solution  is 


s,  =  A  sinhax  +  B  coshax 
h 


(3.39) 


Hence  the  general  solution  is 


s  =  A  sinhax  +  B  coshax  +  — (^?r  ~  wx) 


a 


2  v  2 


(3.40) 


The  boundary  conditions  relating  to  slip  strain 


are 


at  x  =  0 ,  ^  =  0 


(3.41a) 


,  T  ds  _  n 

at  x  =  L ,  =  0 


(3.41b) 
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Differentiating  Equation  3.40  results 
equal  to 

ds  6 

— j —  =  Aacoshax  +  Basinhax  +  (-w) 

dx  2 

a 

By  introducing  the  boundary  condition 
Equation  3.42  becomes 


in  a  slip  strain 

(3.42) 

Equation  3.41a, 


0  =  Aa  -  — ^  w 

a 

from  which 

A  =  -|  w 
a 

Substituting  the  second  boundary  condition, 
into  Equation  3.42  gives 


(3.43a) 


(3.43b) 

Equation  3.41b, 


0 


— 2  coshaL  +  BasinhaL 

a 


(3.44a) 


from  which 


1  -  coshaL 
sinhaL 

Substituting  Equations  3.43b  and  3.44b 
3.40  and  3.42  yields 


(3.44b) 

into  Equations 


s 


sinhax  +  (1 


coshaL) 


coshax 

sinhaL 


+ 


,  wL 
(  2 


wx ) 


(3.45) 


and 

s ' 


coshax  +  (1 


coshaL ) 


s  i  n  h  ax 
sinhaL 


(3.46) 


The  deflection  may  now  be  determined  by 
integration  of  Equation  3.46  using  Equation  3.19(a). 
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Combining  these  two  equations  yields 


-  v"  = 


A  Y 

w  / T  2 .  c  c  w3  i  u 

2E j  (Lx  -  x  )  +  — —  | coshax 

t  t  a 


[■ 


u  t  \  sinhax 

+  (1  -  coshaL)  — : — — 7-  -  1 

sinhaL 


] 


Integrating  twice  yields 

3 


-  v  = 


w  Lx 


2EI 


t 


6 


x 

12 


+ 


A  Y  gw 
c  c 

rt“4 


in  which 


2  2 

+  psinhax  -  a  *  +  D  +  Cx 


p  =  (1  -  coshaL) /sinhaL 


cosh 


(3.47) 


x 


(3.48) 

(3.49) 


condition , 
condition 


0 

from  which 

Subst itut i 


Subjecting  Equation 
v(0)  =  0,  yields  D  = 
v(L)  =  0  yields 


24EI 


T  4  A  Y  gw 

wL  c  c  , 

A  v 


t 


Xta 


C  =  - 


4  3 
a  L 


+ 


a 


24EA  Y  3  '  2 

c  c 


ng  into  Equation  3.48 


3.48  to  the  boundary 
-1.  Imposing  the  boundary 


2.2 

2~~  +  CL) 


(3.50a) 


(3.50b) 


and  rearranging  yields 


-  v 


w 

2EIt 


+ 


A  Y  gw 
c  c 


4“ 


coshax 


4  \  T  3 

x  \  wL  x 

12  ) “  24EI t 

2 

06 

+  psinhax  +  — ^  (Lx 


x2)  -1 


(3.51) 


Equations  3.45,  3.46  and  3.51  give  the 
solution  for  slip,  slip-strain  and  deflection  of  a  simply 
supported  uniformly  loaded  beam.  A  solution  of  the 


differential  equations  for  a  continuous  beam  with  two 
unequal  spans  subjected  to  two  equal  midspan  concentrated 
loads  is  derived  in  Appendix  A.  The  evaluation  of  the 
constants  of  integration  for  this  case  is  so  complex 
that  they  can  only  be  obtained  by  solving  a  matrix 
equation  in  the  computer.  Therefore  solutions  can  only 
be  obtained  once  numerical  values  are  assigned  to  all 
variables . 

3 . 5  Formulation  for  Nonlinear  Load-Slip  Relationships 

The  solutions  in  Section  3.4  and  Appendix  A 
assume  a  linear  relationship  between  shear  flow  and  slip. 
This  assumption  was  introduced  in  Equation  3.11,  where 
it  was  assumed  that 

q  =  Ks  (3.52) 

When  the  slip  is  nonlinear,  Equation  3.52  may  be 
replaced  by  the  general  expression 

q  =  f ( s )  (3.53) 

in  which  f(s)  is  any  nonlinear  function  of  s. 

Yam  and  Chapman  ( 31)  introduced  an 
exponential  relationship  in  the  form 

Q  =  a  (1  -  e-bs )  (3.54) 

to  represent  a  nonlinear  load-slip  relationship  between 
the  total  shear  connector  force,  Q,  and  the  slip,  s. 
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This  may  be  reduced  to  a  shear  flow  relationship  by 
dividing  by  the  connector  spacing,  £.  Then  Equation 
3.52  becomes 

q  =  ^  =  £■  ( 1  ~  e  )  (3.55) 

which  is  a  particular  form  of  the  more  general  relationship 
represented  by  Equation  3.54. 

If  the  relationship  of  Equation  3.53  is  used 
in  place  of  Equation  3.52,  Equation  3.11  becomes 


dP 

dx 


-  f(s)  =  -  q 


(3.56) 


Equation  3.17  becomes 


EA  £  '  -  EA  Y  4)'  =  -  f(s) 

s  0  s  sT 


(3.57) 


and  the  governing  equation  for  slip,  Equation  3.23 


becomes 


s"  -  a  f(s)  =  -  BV 


(3.58) 


in  which,  from  Equation  3.21, 


and 


a 


K 


3  = 


A  A  Y  Y 
EA  _§L  +  ^_s_^ 

A  I 

C  t  t 


A  Y 
s  s 


El ,  A 
t  c 


A 

A 


+ 


A  Y  Y 
s  s  c 


t 


(3.59) 


(3.60) 


Since  Equation  3.58  is  nonlinear,  there  are  no  closed 
form  solutions.  Therefore,  numerical  solution  techniques 
must  be  applied  as  discussed  in  Chapter  IV. 
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FIGURE  3.1 


COMPOSITE  SYSTEM 
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(a) COORDINATE  SYSTEM  &  STRAIN  DISTRIBUTION 


y 


FIGURE  3.2  DEFORMATION  AND  DISPLACEMENT  OF  COMPOSITE  BEAM 
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ii 


FIGURE  3.3  STEP  FUNCTION 
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FIGURE  3.5  CONDITIONS  AT  INFLECTION  POINT 
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CHAPTER  IV 


A  NUMERICAL  SOLUTION  TECHNIQUE 
FOR  COMPOSITE  BEAMS 


4 . 1  Introduction 

Classical  solutions  for  deflections  of 
composite  beams  are  limited  to  simple  sections  with  no 
material  nonlinearities.  Furthermore,  as  may  be  seen  in 
Appendix  A,  the  evaluation  of  constants  of  integration 
becomes  tedious  for  continuous  beams.  Numerical  techni¬ 
ques  offer  an  alternative  in  these  cases. 

A  technique  for  obtaining  numerical  solutions 
for  continuous  composite  beams  with  a  nonlinear  load- 
slip  relationship  is  developed  in  this  Chapter.  The 
governing  equation  for  slip,  Equation  3.23,  is  the  basis 
of  the  numerical  solution  technique.  Beginning  with  an 
initial  assumed  slip,  the  numerical  technique  involves  a 
process  of  iteration  which  is  continued  until  the  slip 
has  an  acceptable  degree  of  accuracy. 
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4 . 2  Numerical  Integration  of  Slip  Equation 


A  three-span  continuous  composite  beam,  shown 
in  Fig.  4.1,  is  divided  into  a  number  of  segments  of 
length  h.  Each  end  of  a  segment  length  is  defined  by 
a  nodal  point,  as  illustrated  in  Fig.  4.1(b).  The 
continuous  beam  is  changed  into  a  simple  beam  primary 
structure  subjected  to  two  sets  of  loadings.  The  first 
loading  consists  of  the  external  loading  while  the  second 
consists  of  the  redundant  reactions.  This  scheme  is 
illustrated  in  Figs.  4.1(c)  and  4.1(d),  respectively. 

As  derived  in  Section  3.5,  the  governing 
second  order  equation  for  slip  is 

s"  -  a2f(s)  =  -  3  V  (4.1) 

Setting  s'  =  v  (which  is  not  to  be  confused  with  v 
of  Chapter  III),  the  second  order  Equation  4.1  may  be 
expressed  as  the  following  two  first  order  simultaneous 
equations : 

s'  =  v  (4.2) 

v'  =  a2f(s)  -  3  V  (4.3) 


Equations  4.2  and  4.3  may  be  solved  by  using  a  Runge- 
Kutta  method  of  integration,  and  iterating  each  equation 
in  turn.  In  applying  the  Runge-Kutta  method  (7) 
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Equations  4.2  and  4.3  are  expressed  as 


=  H(  s  , 

X  ,  v) 

(4.4) 

=  G(s, 

X ,  v) 

(4.5) 

in  which  H  and  G  represent  the  functional  forms  in 
Equations  4.2  and  4.3.  The  recursive  equations  for  the 
fourth  order  Runge-Kutta  solution  of  Equations  4.4  and 
4 . 5  are 


i  +  1  Si  + 

h 

6 

<kls+ 

2k2s+ 

2k3s+ 

k  ) 

4s 

(4.6) 

i+1  Vi  + 

h 

6 

(klv+ 

2k2v  + 

2k3v+ 

k.  ) 

4v 

(4.7) 

Equations  4.6  and  4.7  are  used  alternately  and  the  k. 


and  k.  values  are  determined  as  follows: 

IV 


is 


k.,  =  H(  x  .  ,  s  .  ,  v  .  ) 

Is  v  1  ’  1  ’  1 


k.  =  G(x .  ,  s .  ,  v .  ) 
1  v  v  1  ’  1  ’  1 


(4.8) 

(4.9) 


hk.,  hk1 

k2s  “  H(xi  +  5-  si  +  -2^'  vi  +  (4'10) 


,  hk-,  hk.. 

.  ,  h  ,  Is  ,  lv 

k2v  =  G(xi  +  2’  Si  +  2~ ’  Vi  +  2~ 

v  +  h  +  ^2s  ,  hk2  v 

k3s  H(xi  2’  Si  2  ’  vi  2 


)  (4.11) 

)  (4.12) 


hk. 


,  „ ,  ,  h  ,  ““2s 

k3v  =  G(xi  +  2’  Si  + 


hk 


o  »  v  • 
2  l 


2  v 


+  — ^)  (4.13) 


k4v  =  H(xi  +  h’  si  +  hk3s ’  Vi  +  hk3v}  (4'14) 


k.  =  G(x.  +  h,  s.  +  hk„  ,  v.  +  hkQ  )  (4.15) 

4v  l  l  3s’  l  3v 


' 
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The  initial  values  of  x.  ,  s.  ,  v.  determine  k-i  and  k-, 

from  Equations  4.8  and  4.9  respectively.  The  remaining 

values,  k.  and  k.  ,  where  j  =  2,  3,  4,  are  determined 
js  jv’ 

from  Equations  4.10  through  4.15  successively. 

Equations  4.6  and  4.7  determine  values  of  slip 

and  slip  strain,  respectively.  The  computation  of  slip 

starts  from  the  left  end  of  the  beam  at  x .  =  0  and 

i 

progresses  through  all  nodal  points  to  the  right  end  of 
the  beam  as  shown  in  Fig.  4.1(b).  The  shear  flow  at  any 
nodal  point  in  the  beam  may  be  computed  by  using  either 
a  linear  or  a  nonlinear  relation  for  slip  in  which,  from 
Equation  3.52, 

q  =  Ks  (4.16a) 

or,  from  Equation  3.55, 

q  =  -^-(1  -  e~bs)  (4.16b) 

The  axial  force  in  the  concrete  at  the  right  end  of  the 

beam  is  computed  by  using  the  relationship 

L 

F  =  /  qdx  (4.17) 

0 

If  the  axial  force  computed  at  the  right  end  is  zero,  the 
assumed  initial  slip  value  is  correct  and  hence  the  total 
solution  is  correct.  If  the  axial  force  is  not  zero  at 
the  right  end,  a  correction  is  applied  to  satisfy  the 
boundary  condition 


F(L)  =  0 


(4.18) 
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4 . 3  Correction  and  Iteration  for  Slip 

An  unbalanced  axial  force  in  the  slab  at  the 
right  end  of  the  beam,  as  shown  in  Fig.  4.2a,  creates  an 
unbalanced  moment.  An  equal  and  opposite  balancing 
moment  produces  a  correction  to  the  initial  assumed  slip 
value.  The  general  solution  to  Equation  3.37  for  this 
loading  case  is 


s  =  A  sinhax  +  B  coshax  + 


a 


6  F<Yc  -  V 
L 


(4.19) 


Evaluating  s' ,  and  equating  it  to  zero  for  x  =  0 
requires  that  A  is  zero.  Fig.  4.2(b)  illustrates  the 
moment  correction  to  remove  the  unbalanced  force.  Using 
Equation  3.6,  the  strains  at  the  centroids  of  the  slab 
and  steel  section  at  the  end  of  the  span  are 


£c  =  £0  "  *Yc  +  SL ' 


-F 

EA 


c 


F 


e  =  £_  -  <|>Y  ==  A 

s  O  s  EA 


(4.20) 

(4.21) 


Eliminating  from  Equations  4.20  and  4.21,  the  slip 
strain,  in  terms  of  unbalanced  force  and  curvature,  is 


sl  “  f  (f  +  tt' >  +  *(Yc  -  V 


(4.22) 


From  Equation  3.16  the  slip  strain  can  also  be  expressed 
in  terms  of  curvature  and  the  applied  forces  required 
to  balance  the  unbalanced  moment  as 


sJAY  =  (Y  -  Y  )  +  <t,l. 
Lee  Evc  s  yt 


(4.23) 
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by  eliminating  curvature  from  Equations  4.22  and  4.23 

2- 


si  =  -  = 


L 


A.  I.  +  A  A  (Y  -  Y  ) 
tt  se  e  s 


A  A 
s  c 


I1 


-A  Y  (Y  -Y  )' 
t  c  c  c  s' 


=  X 


(4.24) 


Equation  4.24  defines  the  slip  strain  boundary  condition 
at  the  right  end  of  the  beam  for  the  balancing  moment. 

By  differentiating  Equation  4.19  the  slip  strain  at  the 
right  end  of  the  beam,  where  x  =  L,  is 


s^  =  aB  sinhaL 

By  equating  Equations  4.24  and  4.25, 
evaluated  as 

B  =  A 

asinhaL 


(4.25) 

the  constant  B  is 

(4.26) 


As  the  A  and  B  values  are  now  known,  the  slip  at  the  left 
end,  from  Equation  4.17,  is 

s=B+-^-F(Y  -Y)  (4.27) 

2  v  c  s 
a 

Equation  4.27  represents  the  correction  to  be  applied  to 
the  initial  assumed  slip  at  the  left  end  of  the  beam. 

The  procedure  is  repeated  with  the  corrected  slip  value 
until  the  force  boundary  condition  at  the  right  end, 
Equation  4.20,  is  satisfied.  Values  of  slip  and  slip 
strain  are  then  directly  available  from  Equations  4.6  and 
4.7;  shear  flow  may  be  computed  from  Equations  4.16a  or 
4.16b,  and  axial  force  may  be  computed  from  a  numerical 
integration  of  Equation  4.17. 
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4 . 4  Evaluation  of  Slip  Strain 


Although  the  slip  strain  may  be  obtained  from 
the  Runge-Kutta  solution  as  the  variable  v_^  ,  evaluated  in 
Equation  4.7,  the  slip  strain  used  to  determine  the  shear 
flow  in  the  program  developed  herein  was  obtained  as  the 
derivative  of  a  quadratic  polynomial  passing  through 
three  consecutive  points  on  the  slip  curve.  For  equally 
spaced  points  this  reduces  to  the  standard  second  order 
finite  difference  approximation. 

The  Lagrangian  interpolating  polynomial  f(x), 
passing  through  the  values  f^,  f ^  and  f^  at  three  unequally 
spaced  points  as  illustrated  in  Fig.  4.3(a),  may  be 
expressed  in  terms  of  the  local  x  coordinate,  measured 
from  point  1,  as 

(x-£  x(x-£  ) 

f(x)  -  nr - fi +  f 2  +  und"  f3  <4-28) 


r 


'1  2 


Differentiating  this  expression  yields 


df  =  2(x  ^  Z2  +  £-2x  f  +  ^ 
dx  £_^£  1  £^^2  2  £  £2 


(4.29) 


Evaluating  this  derivative  at  the  three  nodal  points 


gives 

df  ~ ( 2  £^+  £2 ) 


dx 


£-,  +  £9 

f  +  — — —  f 
£x(  £1+£2)  1  £,  ir,  2 


£ 


1 


'1  2 


2  (  £-^+  £2  )  3 


x=0 


(4.30a) 


df 


■£, 


dx  £^(  £-^+  £g )  1 


f„  - 


(£,,-£-,)  £1 

f2  +  £2(V£2)  f3  <4'30b> 
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o 

df  _  z 

dx  1 

x=H 


(Z1+Z2)  Z.+2Z 

V2  f2  +  VW  ^ 


(4.30c) 


These  equations  are  used  to  evaluate  s '  ^ .  In  evaluating 
s1 .  at  node  i,  and  letting  ^  =  h,  Equation  4.30 

becomes 


s'.  = 

l 


l 

2h 


(s  . 


i  +  1 


-  si_i) 


(4.31) 


which  is  the  finite  difference  form. 

At  an  inflection  point  the  section  properties 
may  be  considered  to  be  discontinuous,  as  discussed  in 
Section  3.2.  Let  us  consider  that  the  bending  moment 

diagram  is  as  shown  in  Figure  4.3b  and  the  discontinuity 

in  section  properties  occurs  at  point  0  which  is  located 
at  a  distance  'a'  from  a  node  j  and  at  a  distance  ' b ’ 

from  node  i.  In  terms  of  the  slip  values  at  nodes  j-1, 

j  and  j+1  shown  in  Figure  4.3c,  the  slip  strain  at  points 
j  and  j+1  may  be  computed  by  evaluating  Equations  4.30b 
and  4.30c,  which  become 


and 

s 


j  h( a+h )  Sj-1+  ah  Sj  +  a(a+h)  Sj+1 

a  _  ( a+h )  h+2a 

j+1  h(a+h)  Sj-1  ah  Sk  a(a+h)  Sj+1 


(4.32) 


(4.33) 


Based  on  the  value  of  s’  .  ,  -  ,  the  value  of  s’  .  (Figure 

4.3c)  may  be  evaluated  from  Equation  3.36.  The  standard 

Runge-Kutta  formulae,  Equations  4.6  and  4.7,  may  then  be 

used  to  evaluate  s.  and  s .  ,  .  The  particular  values  of 

l  l+l 


• 
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s’^_^  and  s' . ,  obtained  from  Equations  4.30a  and  4.30b, 
are 

- ( 2b+h )  „  b+h  „  b 

S  i-1  b ( b+h )  Si-1  bh  Si  h(b+h)  si  +  l  (4.34) 


and 

s  '  . 

l 


-h  h-b  b 

b(b+h)  Si-1  h  b Si  h(b+h)  si+l 


(4.35) 


Therefore,  the  slip  strain  can  be  obtained  at  nodes 
adjacent  to  an  inflection  point  by  interpolating  through 
the  values  of  slip  at  these  points. 


4.5  Evaluation  of  Deflections 


Once  the  slip  has  been  determined,  the 
curvature  at  any  point  along  the  length  of  the  beam  may 
be  determined  by  using  Equation  3.16,  since  all  the  terms 
are  known.  The  curvature  is  expressed  by  Equation  3.16 
as 

i<  A  Y 

v"  =  yy  +  -j-£  s'  (4.36) 

Eit  zt 


Defining  v'  =  z  the  above  second  order  equation  may  be 
expressed  as  two  first  order  simultaneous  equations.  Note 
that  v  is  the  displacement  defined  in  Chapter  III. 


v 


z '  = 


M 


El 


+ 


A  Y 
c  c 


s ' 


(4.37) 

(4.38) 


t 


By  employing  a  second  order  Runge-Kutta  method  of 
integration  to  each  equation  in  turn,  Equations  4.37 
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and  4.38  may  be  written  symbolically  as 


v  '  =  G(  z  ,  v ) 

(4.39) 

z '  =  H(z ,  v) 

(4.40) 

from  which 

vi+l  =  1  (klv  +  k2v}  +  vi 

(4.41) 

zi+l  =  1  (klz  +  k2z)  +  Z1 

(4.42) 

Values  of  k.  ,  kn  and  k0  are  evaluated  as 

IV  z  v  z  z 

k  i  =  z . 

lv  1 

M.  AY 

\r  —  1  4.  C  C  ^  1 

1Z  EIt  h  1 

(4.43) 

k0  =  z .  +  k-  x  h 

2v  l  lz 

(4.45) 

M .  ,  A  Y 

k2z  EIt  It  i+l 

(4.46) 

By  substituting  values  from  Equation  4.43  through  4.46  in 
Equations  4.41  and  4.42,  the  deflection  at  every  nodal 
point  is  computed.  Since  the  value  of  v'(0)  is  unknown 
'a  priori' ,  a  value  of  zero  is  first  assumed  before  the 
integration  is  carried  out.  The  deflections  so  obtained 
are  shown  in  Figure  4.4.  To  obtain  the  true  deflection 
at  every  nodal  point,  a  linear  correction  must  be  applied. 

As  illustrated  in  Figure  4.4  the  slope  of 
the  line  joining  A  and  N  is 


0 


L 


(4.47) 
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Therefore,  the  corrected  deflection  at  any  point  is 

A .  =  0 x .  -  6 .  (4.48) 

Thus  a  complete  set  of  deflections  is  obtained. 

The  procedure  described  above  is  valid  for  a 
simply  supported  or  continuous  beam.  However,  to  obtain 
values  for  a  continuous  beam,  it  is  necessary  to  solve  for 
the  indeterminate  reactions,  as  described  in  Section  4.6. 


4.6  Solution  Procedure  for  Continuous  Beams 


The  redundancy  of  a  continuous  beam  may  be 
solved  by  using  the  flexibility  method.  A  three  span 
continuous  beam  with  a  uniform  load  w,  shown  in  Figure 
4.5a,  is  statically  indeterminate  to  the  second  degree. 
The  primary  system  has  been  selected  as  a  simple  beam 
which  deflects  under  applied  loads  as  shown  in  Figs. 

4.5b,  c  and  d.  With  the  redundants  and  deformations  as 
defined  in  these  figures,  the  continuity  equations  can  be 
expressed  as 

5BBRB  +  6BCRC  ~  AB  (4.49) 

6CBRB  +  6CCRC  =  AC  (4.50) 


where  AD  and  A,-,  are  deformations 
5DD)  6^  are  deformations  due  to 

DD  V_,B 

are  deformations  due  to  unit 
are  the  reactions. 


due  to  external  loading; 

unit  load  at  B  and  5  , 

Cd 

load  at  C .  Rn  and  R„ 

Jd  L 
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Equations  4.49  and  4.50  are  expressed  in  the 
matrix  form  as 


PQ 

<o 

1 _ 

■  6bc 

!Rb! 

_  |ab  ) 

(4.51) 

5bc 

6CC 

fRci 

fAc  i 

or 


{R}  =  (A) 


where  F  is  the  flexibility  matrix,  { R } 
vector  and  {A}  deformation  vector.  The 
by  inversion  of  the  flexibility  mat rix 


(4.52) 

is  the  redundant 
reactions,  determined 
are 


{R}  =  {F}  1  {A} 


(4.53) 


The  total  stress  resultants  are  obtained  by  superimposing 
the  stress  resultants  due  to  external  loading  and  those 
due  to  interior  support  reactions. 

The  continuous  beam  of  Fig.  4.5  may  now  be 
solved  by  the  following  procedure.  Assuming  constant 
section  properties,  the  deflections  {A}  for  the  uniform 
load  acting  on  the  primary  structure  are  computed  by  the 
procedure  described  in  Sections  4.2,  4.3  and  4.5, 
Similarly  the  deflections  6T,0  and  6  are  computed  for 
a  unit  redundant  applied  at  B,  and  the  deflections  5 

dL 

and  are  computed  for  a  unit  redundant  applied  at  C. 

The  matrix  {F}  of  Equation  4.52  is  then  formed,  and  the 
redundants  determined  from  equation  4.53. 

If  the  problem  were  linear  a  superposition 


obtained  from 


X 


X0  +  XBRB  +  XCRC 


(4.54) 


in  which  X  represents  any  quantity  (moment,  shear, 

deflection,  slip,  slip  strain,  shear  flow,  etc.),  and 

Xn  ,  X_  and  X_  represent  the  values  obtained  for  that 
U  h>  L 

quantity  from  the  three  basic  solutions  described  in  the 
preceding  paragraphs.  The  solution  is  then  complete  for 
classical  solutions,  such  as  that  of  Plum  and  Horne  (19), 
described  in  Appendix  A. 

Two  types  of  nonlinearity  have,  however,  been 
included  herein.  The  first  is  that  arising  from  the 
cracking  of  the  concrete  in  the  negative  moment  region, 
and  the  second  is  a  nonlinear  slip/shear  flow  relationship. 
In  the  following  paragraphs,  the  continuous  beam  technique 
will  be  considered,  for  the  example  problem  of  Fig.  4.5, 
treating  each  of  the  cracking  effects  separately  and  then 
combining  them  with  the  nonlinear  load-slip  relationship. 

For  the  non  linear it ies  arising  from  the  crackin 
of  the  slab,  the  following  technique  may  be  used.  Super¬ 
position  of  moment  according  to  Equation  4.54  will  produce 
regions  of  negative  moment.  The  locations  of  the  points  of 
inflection  can  be  determined  from  a  linear  interpolation 
of  the  superimposed  moments  between  nodal  points,  resulting 
in  a  number  of  points  along  the  beam  (for  the  present 
example,  four  points)  at  which  a  condition  similar  to  that 
in  Fig.  4.3  arises  For  the  negative  moment  regions  the 
centroid  and  section  properties  of  the  transformed  section 
are  those  of  the  steel  section  and  reinforcing  bars  only. 
Hence  the  discontinuities  indicated  in  Fig.  3.4  arise. 
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The  deflections  {A}  ^  required  to  set  up 
Equation  4.52  are  now  computed  with  these  nonuniform 
section  properties  and  by  considering  the  stress  resultants 
obtained  from  Equation  4.54  with  the  current  values  of  R 

h> 

and  R  .  The  lack  of  compatibility  in  deflection  at  the 
redundants  are  now  obtained  as 

{6A}1  =  {A)x  -  {A}  (4.55) 

Now  the  change  in  redundants  to  eliminate  this  lack  of 
compatibility  is  computed  as 

fJ  {AR)2  =  (6A}1  (4.56a) 

The  improved  estimate  of  reaction  is  then 


{R}2  - 

{  H} 

l  +  CR}2 

(4.56b) 

The 

total  superposition, 

expressed 

by  Equation  4. 

54 ,  is 

now 

repeated  for  the 

new 

react  ions 

unt i 1  ,  for  the 

i  t  h 

iteration 

{AR} 

i 

=  (0) 

(4.57a) 

or 

{ 6  A  } 

i 

=  (0} 

(4.57b) 

When  conditions  4.57a 

or 

4.57b  are 

satisfied  the 

process 

is 

considered  to  have 

converged  an 

d  t he  ' correct ' 

solut i on 

has  been  obtained. 

Since  the  above  iterative  process  depends 
only  on  the  consistency  of  the  final  beam  displacements 
with  the  interior  support  conditions,  it  may  be  applied  to 
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beams  with  a  nonlinear  load-slip  relationship.  The 
primary  difference  is  that  the  iteration  required  to 
determine  the  proper  shear  flow,  described  in  Sect.  4.3, 
which  must  be  carried  for  each  of  the  separate  loading 
conditions,  may  take  longer  to  converge  because  the 
estimated  slip  correction  described  by  Equation  4.27  is 
not  as  good  an  estimate  as  for  the  linear  case. 

4 . 7  Special  Problems  in  Solving  Continuous  Beams 

The  general  technique  developed  herein  for 
the  solution  of  continuous  beams  has  been  described  in 
Sect.  4.6.  In  that  section  it  was  pointed  out  that  at 
an  inflection  point,  there  is  an  abrupt  change  in  section 
properties,  if  it  is  assumed  that  concrete  cannot  resist 
tension.  This  gives  rise  to  the  discontinuities 
illustrated  in  Fig.  3.4.  There  is,  therefore,  a 
discontinuity  in  slip  strain  as  derived  in  Sect.  3.2 
and  expressed  by  Equation  3.36.  It  is  necessary  to 
consider  this  discontinuity  in  the  integration  of  Equations 
4.4  and  4.5,  described  in  Sect.  4.2,  when  solving  continuous 
beam  problems.  This  results  in  the  following  modification 
of  the  procedure  for  the  computation  of  all  of  the  deflections 
required  for  Equation  4.55. 

The  inflection  points  are  located  as  described 
in  Sect.  4.6.  Equation  4.1  is  integrated  by  the  application 
of  Equations  4.2  to  4.15  up  to  node  point  j  of  Fig.  4.3. 
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The  value  of  s  and  v  are  then  determined  for  node  point 
j+1  by  Equation  4.2,  s^  of  Equation  3.36  is  known  and 
hence  s'  is  computed  from  Equation  3.36.  Since  slip  is 
continuous,  the  values  of  s  and  s'  at  mode  i-1  of  Fig.  4.3 
are  now  known,  and  the  values  at  node  i  are  determined  by 
using  ' b '  of  Fig.  4.3  in  place  of  h  in  Equations  4.6  to 
4.15.  The  integration  procedure  now  proceeds  normally  to 
the  next  inflection  point,  when  the  same  procedure  is  used. 
Thus  the  discontinuity  in  slip  strain  is  accounted  for  in 
the  integration  procedure  for  continuous  beams. 


4 . 8  Computer  Program 


A  number  of  computer  programs  v/as  developed 
during  this  investigation.  The  main  program,  developed 
on  the  basis  of  the  numerical  solution  technique,  is 
included  in  Appendix  F.  This  appendix  includes  the  flow 
chart,  description  and  listing  for  the  program,  and  a 
test  problem. 


The  program  is  designed  to  analyse  composite 
beams  of  up  to  six  spans  and  can  handle  any  combination  of 
concentrated  loads  combined  with  a  uniformly  distributed 
loading.  Calculations  may  be  performed  on  the  basis  of: 


( a) 

(b) 


Constant  flexural  stiffness  (EI^), 

Variable  flexural  stiffness  (El  ), 

\  v 


or 

i  .  e  .  , 


different  stiffnesses  in  positive  and  negative  moment 
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regions ,  or 

(c)  Average  flexural  stiffness  (El  ),  i.e., 

v  av 

the  average  of  positive  and  negative  moment  region 
flexural  stiffnesses. 

Computations  for  slip,  slip  strain,  shear 
flow,  force  and  deflections  are  based  on  the  procedures 
described  in  this  Chapter. 


(a)  BEAM  DIVIDED  INTO  SEGMENTS 

y 


A 


(c)  PRIMARY  STRUCTURE 


D 


t  ' 

>  ; 

f  s 

A 

D 

(d)  REDUNDANTS 


FIGURE  4.1  BASIS  FOR  NUMERICAL  PROCEDURE 
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(b)  APPLIED  BALANCE  FORCE  F 


FIGURE  4.2  AXIAL  FORCE  AT  THE  END  OF  THE  BEAM 
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(a)  QUADRATIC  LAGRANGI AN  INTERPOLATION  FUNCTION 


FIGURE  4.3  DISTRIBUTION  OF  SLIP  STRAIN  ACROSS  THE  POINT 

OF  INFLECTION 


FIGURE  4.4  DEFLECTION  AND  ITS  CORRECTION 
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ABC 

(a)  CONTINUOUS  BEAM 


D 


A 


JX 

(b)  PRIMARY  STRUCTURE 


(c)  INFLUENCE  COEFFICIENT  DUE  TO  LOAD 


(d)  FLEXIBILITY  COEFFICIENT  DUE  TO  UNIT  LOAD  AT  B 


(e)  FLEXIBILITY  COEFFICIENT  DUE  TO  UNIT  LOAD  AT  C 


FIGURE  4.5  ILLUSTRATION  OF  FLEXIBILITY  METHOD 


CHAPTER  V 


VERIFICATION  OF  NUMERICAL  TECHNIQUE 
5 . 1  Int  roduct ion 

In  Chapter  IV  a  numerical  technique  was 
established  for  the  solution  of  continuous  composite  beam 
problems  under  the  conditions  of:  (a)  linear  response, 
and  (b)  nonlinear  response  in  which  the  nonlinearity  is 
confined  to  cracking  of  the  concrete  in  the  negative  moment 
region  and/or  a  nonlinear  shear  connector  load-slip 
relationship.  The  basis  for  a  computer  program  to  solve 
such  problems  was  described. 

It  is  the  purpose  of  this  Chapter  (a)  to 
verify  the  numerical  technique  developed  in  Chapter  IV  for 
the  problems  itemized  above  and  (b)  to  establish  the 
adequacy  of  the  model  for  the  computation  of  deflections 
in  simply  supported  and  continuous  composite  beams. 

Since  the  program  does  not  include  the 
effects  of  steel  yielding  or  concrete  crushing,  it  cannot 
be  used  to  compute  deflection  of  composite  beams  at 
ultimate  load.  It  also  does  not  include  the  effect  of 
creep  in  the  concrete.  This  latter  effect,  however,  may 
be  simulated  by  using  a  lower  modulus  of  elasticity  for 
concrete . 
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5 . 2  Program  Verification  for  Linear  Behavior 

5.2.1  Simply  Supported  Beams 

The  classical  solution  for  a  simply  supported 
uniformly  loaded  beam  with  a  flexible,  but  linear  elastic, 
shear  connection  is  derived  in  Sect.  3.4.  A  comparison  of 
the  numerical  results  obtained  from  the  program  described 
in  Chapter  IV  with  this  classical  solution  is  contained  in 
Figs.  5.1  to  5.3.  The  properties  of  the  chosen  cross- 
section  are  contained  in  Column  (1)  of  Table  5.1.  The 
classical  solution  for  a  simply  supported  beam  subjected 
to  a  concentrated  load  is  included  in  Appendix  B  and  the 
results  obtained  from  numerical  analysis  are  compared  with 
the  closed  form  solution. 

It  can  be;  seen  from  Figs.  5.1  to  5.3  that, 
based  on  linear  load-slip  behavior  of  shear  connector, 
the  results  of  the  numerical  solution  for  slip,  slip  strain 
and  deflection  (at  a  load  of  30  kips)  are  indistinguishable 
from  the  values  obtained  in  the  classical  solution. 
Nonlinear  load-slip  behavior  for  the  shear  connectors 
increases  the  magnitude  of  slip  and  deflection  throughout 
the  length  of  the  beam. 

In  addition  to  the  solution  for  the  linear 
elastic  problem,  solutions  for  complete  interaction, 
obtained  by  setting  the  shear  connector  stiffness  to  a 
very  high  value,  and  for  nonlinear  load-slip  behavior 
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are  also  shown  on  the  figures.  In  Fig.  5.3  deflection, 
including  the  effect  of  shearing  deformation,  is  also 
shown  for  comparison  purposes. 

This  example  verifies  that  the  results 
obtained  by  the  numerical  technique  for  simply  supported 
beams  are  consistent  with  the  linear  elastic  classical 
solution . 

5.2.2  Continuous  Beams 

The  ability  of  the  computer  program  (based 
on  numerical  technique)  to  properly  analyze  a  continuous 
beam  with  flexible  shear  connectors  but  linear  elastic 
response  is  verified  in  this  Section  by  solving  the 
problem  which  was  the  subject  of  investigation  by  Plum  and 
Horne  (19).  The  problem  is  illustrated  in  Fig.  5.4.  Plum 
and  Horne  solved  this  problem  by  determining  the  general 
solution  to  the  governing  equations,  derived  in  Chapter 
III,  and  imposing  the  proper  boundary  conditions  in  order 
to  evaluate  the  constants  of  integration.  For  this  solution 
constant  section  properties  throughout  the  length  of  the 
beam,  that  is,  concrete  cracking,  is  ignored  in  the  nega¬ 
tive  moment  region. 

Since  Plum  and  Horne  did  not  provide  details 
of  their  solution,  the  classical  solution  for  a  two-span 
continuous  composite  beam  has  been  derived  in  Appendix  A. 

A  computer  program  to  evaluate  slip,  slip  strain  and 


77 


deflection  consistent  with  the  classical  solution  derived 
in  Appendix  A  was  written,  and  the  problem  of  Fig.  5.4 
was  solved  in  this  manner.  The  properties  of  the  beam  are 
tabulated  in  Column  (2)  of  Table  5.1.  The  solution  for 
slip,  slip  strain,  slab  force,  and  deflection  obtained  by 
the  numerical  technique  of  Chapter  IV  are  compared  with 
the  classical  solution  in  Figs.  5.5  to  5.8,  respectively. 

The  only  significant  discrepancy  between  the 
classical  solution  and  the  numerical  solution  is  that  for 
slip  strain  in  the  region  of  the  concentrated  loads  and  the 
reactions.  This  discrepancy  appears  to  arise  because  the 
numerical  solution  technique  computes  slip  strain  by 
averaging  values  obtained  from  the  two  element  lengths  on 
each  side  of  the  nodal  points.  This  technique  does  not 
permit  the  accurate  evaluation  of  slip  strain  at  the  point 
of  application  of  a  concentrated  load.  Nevertheless,  the 
correspondence  between  slip  curves  indicates  that  the 
distribution  of  shear  flow  (directly  proportional  to  slip 
in  a  linear  problem)  predicts  an  accurate  determination  of 
the  distribution  of  internal  forces.  The  deflections  are 
also  reliable. 


It  is  concluded  from  this  problem  that  the 
numerical  techniques  incorporated  into  the  computer  program 
are  suitable  for  the  solution  of  composite  beam  problems 
involving  linear  load-slip  response  of  shear  connectors 
and  concentrated  loads. 
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5 . 3  Program  Verification  for  Nonlinear  Behavior 

Any  solution  which  includes  nonlinearities 
for  composite  beams  must  be  based  on  a  numerical  technique. 
Few  of  these  solutions  are  available.  However,  Yam  and 
Chapman  (31)  have  published  numerical  solutions  based  on 
predictor  corrector  technique  for  continuous  composite  beams 
which  include  the  effects  of  steel  yielding,  concrete 
tensile  cracking,  concrete  crushing  and  a  nonlinear  expo¬ 
nential  relationship  between  force  and  slip  in  the  shear 
connectors . 

The  primary  objective  of  this  thesis  is  to 
investigate  deflections  of  composite  beams  under  service 
loads.  Nonlinearities  due  to  yielding  of  the  steel  and 
crushing  of  the  concrete  have  no  effect  at  this  level  of 
load.  Therefore,  the  simpler  technique  developed  herein 
is  appropriate  for  the  present  study.  In  order  to  demon¬ 
strate  this,  the  solution  obtained  from  the  present 
numerical  technique  is  compared,  at  service  load,  with  the 
full  nonlinear  solution  of  Yam  and  Chapman. 

Yam  and  Chapman's  test  beam  is  shown  in 
Fig.  5.9  and  its  properties  are  shown  in  Column  (3)  of 
Table  5.1.  Their  solution  for  slip  and  deflection,  at  a 
service  load  of  7.4  tons,  is  shown  in  Figs.  5.10  and  5.11, 
respectively,  where  they  are  compared  with  the  present 


numerical  solution. 
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The  comparison  of  slip  in  Fig.  5.10  indicates 
that,  although  there  are  some  differences  in  the  vicinity 
of  the  concentrated  loads,  the  characteristics  of  the  curves 
are  very  similar.  In  particular,  the  maximum  slip  values 
are  essentially  the  same  and  occur  at  approximately  the  same 
locations.  The  comparison  of  the  service  load  deflections 
for  these  two  solutions  is  shown  in  Fig.  5.7  where  the 
values  are  seen  to  be  identical  for  all  practical  purposes. 

The  increase  in  deflection  produced  by  the  flexible  shear 
connectors  and  concrete  cracking  may  be  seen  by  comparing 
values  with  the  complete  interaction  deflections,  shown  also 
in  Fig.  5.11  and  which  are  taken  from  the  paper  by  Yam  and 
Chapman.  From  Fig.  5.11  it  is  evident  that  shear  deformation 
has  not  been  included  in  Yam  and  Chapman's  analysis.  The 
degree  of  nonlinearity  developed  in  the  load-slip  relationship 
for  this  problem  may  be  seen  by  determining  the  position  at 
which  the  maximum  slip  of  Fig.  5.6  is  located  on  a  load-slip 
or  shear-flow-slip  curve.  This  is  shown  in  Fig.  5.12, 
from  which  it  may  be  concluded  that  there  is  some  nonlinearity 
associated  with  the  load-slip  relationship  at  service  load 
for  this  problem. 

Figs.  5.10  to  5.12  appear  to  support  the 
argument  that  an  analysis  in  which  the  non lineari t ies  are 
confined  to  cracking  of  the  concrete  and  the  load-slip 
relationship  is  capable  of  producing  an  adequate  simulation 
of  the  behavior  of  composite  beams  at  service  loads. 
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Furthermore,  this  establish  that  the  program  developed 
herein  has  incorporated  these  nonlinearities  in  such  a  way 
that  the  results  are  consistent  with  a  more  sophisticated 
treatment  of  the  problem.  Thus,  it  is  considered  that  the 
reliability  of  the  computational  technique  is  adequate  to 
draw  conclusions  relative  to  the  effects  on  service  load 
deflections  arising  from  variations  on  stiffness  parameters. 

5 . 4  Comparisons  with  Test  Results 

5.4.1  Comparison  with  Hamada  and  Longworth  (  9  ) 

Hamada  and  Longworth  (  9  ),  in  an  experimental 
and  analytical  investigation  into  the  behavior  of  continuous 
composite  beams,  concluded  that  shear  deformation  may 
contribute  significantly  to  the  deflection  of  composite 
beams.  Load  deflection  computations  for  one  of  their  beams 
(CBI)  are  shown  in  Fig.  5.13.  The  properties  used  in  the 
analysis  of  this  beam  are  summarized  in  Column  (4)  of 
Table  5.1.  The  computation  of  service  load  for  this  beam 
is  given  in  Appendix  C  and  the  results  are  summarized  in 
Column  (2)  of  Table  5.2.  Very  good  agreement  is  obtained 
with  Hamada  and  Longworth ' s  results  up  to  70  kips  using 
nonlinear  load-slip  behavior  of  shear  connectors  and 
including  the  shear  deformations.  It  is  apparent  that, 
for  the  shear  connector  stiffness  used  in  this  analysis,  it 
is  necessary  to  include  shear  deformations  to  obtain 
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reasonable  agreement  with  their  experimental  results. 
Therefore,  shear  deflection  will  be  included  in  many  of 
the  parameter  studies  of  Chapter  VI. 

5.4.2  Comparison  with  Chapman  and  Balakrishnan  (  5  ) 

Chapman  and  Balakrishnan  (  5 )  carried  out  an 
extensive  testing  program  of  simply  supported  composite 
beams  with  various  types  of  shear  connectors.  A  comparison 
with  one  of  their  tests  (Beam  A6)  is  shown  in  Fig.  5.14. 

The  beam  properties  are  summarized  in  Column  (5)  of  Table 
5.1.  The  service  load  is  computed  in  accordance  with  the 
procedure  of  Appendix  C  and  the  computation  is  summarized 
in  Column  (3)  of  Table  5.2. 

In  contrast  to  Hamada  and  Longworth's  beam 
it  should  be  noted  that  better  agreement  is  obtained  with 
these  results  when  linear  load-slip  behavior  is  assumed. 
This  appears  to  contradict  the  conclusion  obtained  from 
the  comparison  in  Sect.  5.4.1.  However,  there  are 
uncertainties  involved  in  evaluating  shear  connector 
stiffness  for  Hamada  and  Longworth's  beam.  Hence  solutions 
with,  and  without,  shear  deformations  will  be  presented 


in  Chapter  VI . 
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5 . 5  Conclusion 

The  numerical  technique  developed  in  Chapter 
IV  predicts  the  behavior  of  composite  beams  with  very  good 
accuracy  up  to  service  loads,  and  hence  the  validity  of 
the  model  is  established. 


TABLE  5.1  PROPERTIES  OF  COMPOSITE  BEAMS 
FOR  TEST  PROBLEMS 


Simply 

Plum  and 

Yam  and 

Chapman  and 

Hamada  and 

Supported 

Horne ’ s 

Chapman ’ s 

Balakrishnan 1 s 

Longworth 1  s 

Beam 

Beam 

Beam 

Beam 

Beam 

Beam  Designation 

BSS12  x  6L 

W12  x  31 

BSB  108 

BSS12  x  6L 

W12  x  31 

Total  Depth=d  in. 

18.000 

16.090 

8.375 

18.000 

16.090 

Concrete  Slab 

be  ( in . ) 

48.000 

48.000 

19.000 

48.000 

48.000 

tc  ( in . ) 

6.000 

4.000 

2.375 

6.000 

4.000 

Ec  (ksi) 

3. 550xl03 

3. 750xl03 

4.002xl03 

3 . 550x 103 

4 . 527xl03 

f'c  (psi) 

3440.0 

4000.0 

6900.0 

3440.0 

5577.0 

Reinforcement 

Agr  (sq. in. ) 

- 

1.571 

0.690 

- 

1.600 

d  ( in.  ) 

- 

3.000 

1 . 750 

2.000 

f  (ksi) 

— 

50.000 

46.600 

- 

59.300 

yr  v 

E  ( kst ) 

30. OOOxlO3 

30. OOOxlO3 

_ 

30. OOOxlO3 

sr 

Steel  Section 

ds  ( in .  ) 

12.000 

12.090 

6.000 

12.000 

12.090 

Ag  (sq. in. ) 

13.000 

9.130 

3.530 

13.000 

9.130 

bf  (in. ) 

6.000 

6.530 

3.000 

6.000 

6.530 

tf  ( in.  ) 

0.717 

0.465 

0.377 

0.717 

0.465 

t  ( in .  ) 
w  v  ' 

f  f  (ksi) 
y  f  ' 

0.400 

0.265 

0.230 

0.400 

0.265 

34.832 

- 

- 

34.832 

40.500 

f  (ksi) 

38.820 

44.000 

46.600 

38.820 

46.500 

y  w 

Es  (ksi) 

31 . 140xl03 

30. OOOxlO3 

30 . 016x 103 

31. 140x  103 

30. 200xl03 

Shear  Connector 

Type 

Headed  Type 

Size 

3/4"x4" 

3/4”x3" 

3/8"x2 ! 

3/4"x4” 

3/4"x3” 

4 

K  (psi) 

150. OOOxlO4 

150. OOOxlO4 

30. 440xl04 
25. 360xl02 

14.400x10 

19 . 250x1(j 

2 

a  (psi) 

18 . 200xl02 

- 

18.200x10“ 

21 . 290x10 

b 

79.160 

- 

120.000 

79.000 

90,418 

Qu  (kips) 

28.000 

28.760 

7.280 

28.716 

28.716 

J 


■ 
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TABLE  5.2  TEST  PROBLEM  COMPUTED  VALUES 


Hamada  and 
Longworth ' s 
Beam 

Chapman  and 
Balakrishnan ' s 

Beam 

Beam  Designation 

W12  x  31 

BSS12  x  6L 

Shear  Connectors 

Type 

Headed  Stud 
(Paired) 

Headed  Stud 
(Paired) 

Size 

3/4"x3" 

3  /  4 ' 1 4  " 

na 

16 

16 

nb 

16 

— 

nc 

8 

— 

Qu  (kips) 

28.716 

28.000 

C  =  .85f'bctc 

910.170 

842.112 

(kips) 

T  =  A  f  (kips) 

387.456 

469.670 

na  Qu  (kips) 

459.456 

448.000 

nb  Qu  (kips) 

459.456 

- 

nc  Qu  (kips) 

229.728 

— 

MBu  (ft--kiPs) 

291.775 

393.000 

MDu  (ft. -kips) 

203.147 

- 

Pu  (kips) 

129.384 

83.463 

P  (kips) 

86.250 

55.650 

w 
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FIGURE  5.1  COMPARISON  OF  CLASSICAL  AND  NUMERICAL  VALUES  FOR  SLIP 
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FIGURE  5.3  COMPARISON  OF  CLASSICAL  AND  NUMERICAL  VALUES  FOR  DEFLECTION 


88 


W  W 


FIGURE  5.4  PLUM  AND  HORNE'S  TEST  BEAM 


CLOSED  FORM  (PLUM  &  HORNE) 
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FIGURE  5.5  COMPARISON  OF  SLIP  VALUES  (PLUM  AND  HORNE  AND  NUMERICAL  SOLUTION) 
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FIGURE  5.6  COMPARISON  OF  SLIPSTRAIN  VALUES  (PLUM  AND  HORNE  AND  NUMERICAL  SOLUTION) 


CLOSED  FORM  (PLUM  &  HORNE) 
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FIGURE  5.7  COMPARISON  OF  AXIAL  FORCE  VALUES  (PLUM  AND  HORNE  AND  NUMERICAL  SOLUTION) 
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FIGURE  5.8  COMPARISON  OF  DEFLECTION  VALUES  (PLUM  AND  HORNE  AND  NUMERICAL  SOLUTION) 
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FIGURE  5.9  YAM  AND  CHAPMAN'S  TEST  BEAM 
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FIGURE  5.10  COMPARISON  OF  SLIP  VALUES  (YAM  AND  CHAPMAN  AND  NUMERICAL  SOLUTION) 
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FIGURE  5.12  NON-LINEAR  SHEAR  FLOW  AND  SLIP  RELATIONSHIP  (YAM  AND  CHAPMAN) 
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FIGURE  5.13  COMPARISON  OF  DEFLECTION  VALUES  (HAMADA  AND  LONGWORTH  AND  NUMERICAL  SOLUTION) 
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CHAPTER  VI 


APPLICATIONS 


6 . 1  Int  roduct ion 

In  Chapter  V  the  capabilities  of  the  numerical 
solution  developed  in  Chapter  IV  were  compared  with  closed 
form  solutions  and  published  results.  It  was  concluded 
that  the  program  developed  herein  is  adequate  to  reliably 
predict  deflections  of  simply  supported  and  continuous 
composite  beams  for  a  variety  of  loading  conditions,  and 
that  the  nonlinearities  included  in  the  program  capture 
the  essential  aspects  of  behavior  at  service  loads.  It 
is  necessary  to  attempt  to  draw  some  general  conclusions 
with  respect  to  the  effect  of  shear  connector  flexibility 
and  nonlinear  behavior  on  service  load  deflections. 

In  this  Chapter  an  evaluation  of  the  effect  of 
various  parameters  on  service  load  deflections  is  attempted. 
For  simply  supported  beams  which  respond  in  the  linear 
elastic  range  at  service  loads,  nondimensional  parameters 
are  derived  from  the  closed  form  solution.  This  permits 
a  quantitative  evaluation  of  increments  in  deflections 
relative  to  those  which  would  occur  in  a  beam  with  complete 
interaction.  Assuming  the  same  nondimensional  parameters 
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control  the  behavior  of  continuous  beams,  the  numerical 
technique  developed  herein  is  used  for  the  evaluation  of 
deflection  in  continuous  beams.  However,  for  nonlinear 
response  no  valid  method  of  nondimensionalizing  the  results 
has  been  found.  Therefore,  it  has  only  been  possible  to 
examine  the  effects  of  the  parameters  on  particular  example 
problems,  with  the  inference  that  the  same  qualitative 
type  of  result  occurs  in  other  beams  of  similar  type. 


6 . 2  Deflections  of  a  Linear  Elastic  Simply  Supported 

Composite  Beam  Subjected  to  Uniform  Load 


The  solution  of  the  differential  equation  for 

a  simply  supported  beam  subjected  to  uniform  load  was 
determined  in  Sect.  3.4.  Evaluating  Equation  3.51  at 


x  =  L/2  yields  the  maximum  deflection  in  the  span  as 

4 

D  W1 

-  V. 


L/2 


r.  j  ...  AY  gw  .  T 

_JL_  wL_  +  c  c /  h  oL 

384  El.  T  4  '  OStl  3 


t 


Zta 


2  2 

,  (1  -  coshaL)  .  ,  aL  ,  a  i 

+  - : — r—= smh  +  — 5 —  -  1 

sinhcil  2  8 


(6.1) 


Considering  downward  deflection  as  positive,  and 
recognizing  that  the  first  term  is  the  deflection  of  a 
beam  with  complete  interaction,  Equation  6.1  may  be 


written  as 


-  —  =  1  + 
60 


384  EAcYcB 
5  a* * * 4L4 


1  - 


2.  2 

a  L 

“ 8““ 


(1  -  coshaL) 


sinhaL 


.  ,  aL  ,  aL 

smh  ~2~  -  cosh 


(6.2) 
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in  which  6^  is  the  midspan  deflection  (excluding  shear 
deformation)  of  the  composite  beam  and  6^  is  the  midspan 
deflection  of  the  same  beam  (due  to  bending  only)  assuming 
full  interaction. 

Recalling  the  definitions  of  a  and  3  from 
Equations  3.21  and  3.22,  Equation  6.2  may  be  written  in 
terms  of  two  nondimensional  ratios  as 


■j—  -  1  + 

60 


384 

c-  r  2 

5£n 


1  -  -  cosh 


a.-  cosh  /p)  sinh  /p; 

smh  /|r T  2 


in  which  the  following  definitions  apply. 

A,  Y  Y 

 t  s  c 


I+.  +  A ,  Y  Y 
t  t  s  c 


and 


n  =  - 


KitL‘ 


EA  A  Y  Y 
s  c  s  c 


(6.3) 


(6.4) 


(6.5) 


Note  that  £  and  p  are  always  positive.  It  can  be  seen 
that  the  ratio  E,  is  dependent  only  on  the  geometry  of  the 
cross-section,  whereas  the  ratio  q  depends  on  the  geometry 
of  the  cross-section,  the  span,  the  shear  connector  stiffness 
factor  (K),  and  the  elastic  moduli  of  steel  and  concrete. 

The  right  hand  side  of  Equation  6.3  represents 
a  factor  by  which  6q  can  be  multiplied  in  order  to  determine 
the  maximum  delfection  of  a  composite  beams  with  any  shear 
connector.  It  is  universally  applicable  to  the  type  of 
span  considered  and  can  form  the  basis  for  a  simple  design 
aid.  To  determine  the  range  of  the  nondimensional  parameters 
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within  which  most  composite  beams  would  fall,  and 
were  computed  for  a  number  of  typical  sections.  By 
arbitrarily  varying  these  factors  over  this  range,  the 
plot  of  Fig.  6.1  was  produced.  This  figure  indicates  that 
shear  connector  flexibility  may  have  a  considerable  effect 
on  composite  beam  deflection.  With  reduction  in  shear 
connector  stiffnes,  that  is,  with  partial  shear  connection, 
the  shear  connector  modulus  K  is  reduced.  Due  to  this 
reduction  in  Equation  6.5,  n  decreases  and  hence  the 
deflection  due  to  slip  increases.  (See  Fig.  6.1.) 
Depending  on  the  properties  of  the  cross  section,  the 
geometry  of  the  beam  and  the  shear  connector  modulus, 
this  deflection  due  to  slip  can  be  as  high  as  45%  of  the 
deflection  due  to  flexure. 


6 . 3  Evaluation  of  Shear  Deflection  for  a  Simply 

Supported  Composite  Beam  Subjected  to  Uniform  Load 


as  (  23) 


The  shear  strain  in  a  beam  may  be  expressed 

(6.6) 


Y  =  v'  , 
s(x) 


K  .  V 
sh 


GA 


t 


in  which  v  is  the  shear  deflection  and 
s 


A 


K 


t 


sh 


/ 


Q‘ 


dA 


(6.7) 


t  A 

in  a  property  of  the  cross-section.  For  a  simply 
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supported  beam,  symmetrically  loaded, 


6 

s 


V 


s(x) 


dx 


V  dx 


(6.8) 


Shear  deflection  is  therefore  simple  to  evaluate  provided 

K  ,  is  known, 
sh 


hand  side 
ratio  of 
load  as 


or 


Recognizing  that  the  integral  on  the 
of  Equation  6.8  is  the  moment  at  midspan 
shear  deflection  to  6^  can  be  written  for 


K  ,  _ 2  384EI , 

sh  wL  t 


0 


GA 


t  8 


5wL 


6s  =  48  EKshrt 
^  5  GAtL2 


right 
,  the 
uniform 

(6.9a) 

(6.9b) 


Unfortunately  this  ratio  cannot  be  expressed  directly  in 

terms  of  s  and  3  and,  therefore,  it  is  probably  best  to 

°s  ^  f 

compute  -r—  separately  and  add  it  to  7 —  if  shear  deflections 

ao  60 

are  to  be  combined  with  flexural  deflections. 

However,  since  the  evaluation  of  K  ,  by 

s  in 

Equation  6.7  is  not  simple,  its  value  has  been  computed 
for  a  number  of  cross-sections  and  is  plotted  against  the 
width  of  the  concrete  slab  in  Fig.  6.2.  K  ,  for  various 
wide  flange  steel  sections  with  different  concrete  slab 
widths  and  for  thickness  t  =  6"  and  4"  are  tabulated  in 
Table  6.1.  It  is  observed  that  Kg^  values  for  the  beams 
lie  between  3.317  and  10.229.  Kg^  for  the  beams  not 
included  in  the  table  may  be  derived  from  the  upper  and 


104 


lower  bound  values  of  the  beam  of  the  same  designated 
rolled  section  using  the  linear  variation  in  depth.  Using 
these  shear  form  factors,  shear  deflection  for  simply 
supported  beams  can  be  computed  with  ease  for  practical 
design  purposes. 

The  magnitude  of  shear  deflection  for  a 
specific  simply  supported  uniformly  loaded  beam  has  been 
shown  in  Fig.  5.3  which  indicates  that  it  is  not  negligible. 

6 . 4  Influence  of  Type  of  Loading,  Shear  Connection  and 
Load-Slip  Characteristics  of  the  Shear  Connectors 
and  Shear  Deformation  on  Deflections 

6.4.1  Simply  Supported  Beam  Subjected  to  Uniformly 
Distributed  Load  and  a  Concentrated  Load 

Beam  U-6  tested  by  Chapman  and  Balakrishnan 
(  5  )  was  modeled  with  the  same  parameters,  and  the  load- 
deflection  curve  up  to  80  kips  (service  load)  was  obtained 
using  the  numerical  technique  discussed  in  Chapter  IV. 

Fig.  6.3  shows  the  comparison  of  deflections  for  full 
shear  connection  and  partial  shear  connection  (50%  of  full 
shear  connection)  considering 

(a)  the  linear  and  nonlinear  load-slip 
characteristics  for  the  shear  connector 

(b)  the  effect  of  shear  deformations. 
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It  can  be  seen  that  each  effect,  namely  slip,  nonlinear 
load-slip  characteristics  and  the  effect  of  shear  defor¬ 
mation  progressively  increase  deflection  and  slip  and  that 
deflection  and  slip  for  beams  with  partial  shear  connection 
are  somewhat  larger  than  for  those  with  full  shear  con¬ 
nection.  The  total  deflection  of  the  composite  beam 
considering  the  effects  of  nonlinear  load-slip  character¬ 
istics,  shear  and  partial  shear  connection  is  approximately 
twice  that  of  the  bending  deflection  of  the  beam  neglecting 
the  above  effects.  A  similar  trend  is  observed  in  the  case 
of  beams  subjected  to  concentrated  load  and  the  deflections 
are  significantly  greater  than  those  for  uniform  load.  The 
percentage  increase  in  deflection  as  well  as  end  slip  from 
full  shear  connection  to  partial  shear  connection  is  given 
in  Table  6.2  and  Table  6.3  respectively.  Refering  to 
Table  6.2,  it  can  be  seen  that  the  maximum  percentage 
increase  in  deflection  is  found  when  both  nonlinear  load- 
slip  characteristics  for  the  connectors  as  well  as  shear 
deformations  are  included.  The  properties  of  the  beams  are 
given  in  Columns  1  and  2  of  Table  6.5. 

Fig.  6.4  shows  the  comparison  of  load  deflection, 
load-slip  curves  for  full  and  partial  shear  connection 
considering  a)  linear  and  nonlinear  load-slip  characteristics 
of  the  shear  connectors  and  b)  the  effect  of  shear  deformation 
for  simply  supported  beam  subjected  to  uniformly  distributed 
load.  Similar  load  deflection  and  load  slip  curves  are 
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drawn  in  Fig.  6.5  and  6.6,  when  the  beam  is  subjected  to 
concentrated  load.  From  the  above  figures,  it  can  be 
observed  that  the  deflections  and  slip  of  the  beam  are 
increased  due  to  nonlinear  load-slip  characteristic  of 
the  shear  connectors  and  the  shear  deformations.  The 
effect  of  partial  shear  connection  on  deflection  and  slip 
is  larger  than  that  of  full  shear  connection. 


6 . 5  Deflection  of  a  Continuous  Beam  Subjected  to 
Uniform  Load 


Since  a  closed  form  solution  for  the  deflection 
of  a  continuous  beam  considering  nonlinear  load-slip 
characteristics  of  shear  connector  is  not  available,  it 
is  not  possible  to  derive  an  expression  of  the  form  of 
Equation  6.3  to  determine  the  effect  of  flexible  shear 
connectors  on  continuous  beam  deflections.  However,  it 
is  reasonable  to  assume  that  the  same  nondimensional 
factors  which  control  the  deflection  os  simply  supported 
beams  also  control  the  deflection  of  continuous  beams. 

A  universally  applicable  plot  similar  to  that  of  Fig.  6.1 
can  then  be  constructed  from  the  results  of  numerical 
solutions . 


Fig. 

with  two  equal 
of  results  are 


6.7  is  such  a  plot  for  a  continuous  beam 
spans  subjected  to  uniform  load.  Two  sets 
shown,  namely,  those  with  zero  stiffness 
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of  the  concrete  in  the  negative  moment  region  (variable  I) 
and  those  with  constant  moment  of  inertia. 

Fig.  6.8  is  identical  to  Fig.  6.7  except  that 
it  includes  the  effect  of  shear  deformation.  Recalling 
the  discussion  from  Sect.  6.3,  it  was  pointed  out  that 
shear  deflection  is  not  directly  a  function  of  £  and  n . 
Nevertheless,  the  smooth  variation  of  the  curves  indicates 
that,  in  the  absence  of  better  information,  Fig.  6.8 
could  give  a  reasonable  indication  of  the  deflection  of 
continuous  beams  including  the  shear  effect.  For  a  more 
accurate  computation,  shear  deflections  can  be  super¬ 
imposed  on  the  results  from  Fig.  6.7,  as  indicated  in 
Sect .  6.3. 

A  comparison  between  Fig.  6.8  and  6.7  indicates 
the  very  significant  increase  in  deflection  due  to  the 
effect  of  shear  deformations  in  continuous  beams.  This 
fact  has  previously  been  pointed  out  by  Hamada  and 
Longworth  (9).  Since  the  shear  deflection  of  a  continuous 
beam  subjected  to  a  concentrated  load  is  approximately  the 
same  as  that  of  a  simple  beam,  the  increase  in  deflection 
ratio  results  primarily  from  the  fact  that  the  flexural 
deflection  in  a  continuous  beam  is  approximately  one- 
quarter  of  that  for  a  simple  beam  and,  hence,  the  shear 
deflection  is  a  higher  fraction  of  the  flexural  deflection 
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0 . 6  Deflection  of  a  Two-Span  Continuous  Beam 
Subjected  to  Concentrated  Load 

The  properties  of  the  beam  are  given  in  Column  3 
of  Table  6.5.  Figs.  6.9  and  6.10  show  the  comparison  of 
load-deflection  and  load-slip  curves  for  full  and  partial 
shear  connection  considering  a)  linear  and  nonlinear 
load-slip  characteristics  for  the  shear  connector  and 
b)  the  effect  of  shear  deformations.  The  figures  are 
based  on  constant  moment  of  inertia.  Similar  curves  are 
drawn  in  Figs.  6.11  and  6.12  for  a  beam  with  variable 
moment  of  inertia.  In  continuous  beams  it  may  also  be 
observed  that  each  effect,  namely  slip,  nonlinear  load- 
slip  behavior  of  the  shear  connectors  and  the  effect  of 
shear  deformation  increases  deflection  and  slip.  Partial 
shear  connection  causes  larger  deflection  than  that 
obtained  for  full  shear  connection.  The  deflection 
comparison  of  continuous  beams  at  service  load  for  partial 
and  full  shear  connection  is  given  in  Table  6.4. 

6 . 7  Influence  of  Negative  Moment  Cracking  in 

Continuous  Beams 

While  the  figures  in  Sect.  6.5  give  the  influence 
of  negative  moment  cracking  on  deflections,  they  give  no 
insight  into  the  effect  of  this  cracking  on  the  distribution 
of  slip  and  slip-strain  along  the  length  of  the  beam. 
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The  effect  of  cracking  on  these  distributions  is  examined 
for  a  typical  beam  in  this  Section. 

The  example  beam  has  the  properties  listed  in 
Column  4  of  Table  6.5.  The  distribution  of  deflection, 
slip,  slip-strain  and  force  in  the  slab  are  shown  in 
Figs.  6.13  to  6.16,  respectively,  for  the  beam  with  and 
without  cracking.  It  can  be  seen  that  the  effect  of 
cracking  is  to  increase  the  deflections  (Fig.  6.13), 
increase  slip  in  the  positive  moment  region  (Fig.  6.14), 
introduce  discontinuities  in  the  slip-strain  relationship 
(Fig.  6.15),  and  increase  the  slab  force  in  the  positive 
moment  region  (Fig.  6.16). 

6 . 8  Analyses  of  a  Three-Span  Beam  for  Shored  and 

Unshored  Construction 

The  example  problem  is  a  beam  with  three  equal 
spans  subjected  to  uniform  load  with  properties  listed 
in  Column  5  of  Table  6.5.  The  live  load  to  dead  load 
ratio  is  2.25.  In  the  analysis  of  unshored  construction 
all  the  dead  load  is  assumed  to  be  carried  by  the  steel 
beam  and  the  live  load  by  the  composite  action.  Full  shear 
connection  and  linear  load-slip  behavior  is  assumed. 
Deflection,  slip,  slip  strain  and  force  distribution  in 
Fig.  6.17  to  6.20.  The  effect  of  unshored  construction 
is  to  increase  the  deflection,  decrease  slip  in  both 
positive  and  negative  moment  regions  and  decrease  slip 
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strain  as  well  as  introduce  discontinuities  in  the  slip- 
strain  relationship  and  decrease  the  slab  force  in  both 
positive  and  negative  moment  regions. 

6 . 9  Conclusions 

In  this  Chapter  an  evaluation  of  the  effect  of 
various  parameters  such  as 

(a)  partial  and  full  shear  connection 

(b)  linear  and  nonlinear  load  slip  behavior 

(c)  shear  deformation 

(d)  constant  and  variable  moment  of  inertia,  and 

(e)  shored  and  unshored  construction 

on  service  load  deflections  have  been  studied.  The 
computer  program  reliably  predicts  deflection  and  slip 
for  simply  supported  and  continuous  beams  under  service 
loads.  The  shear  form  factors  for  various  sections  with 
different  widths  and  thicknesses  of  concrete  slab  are 
tabulated  for  practical  design.  Shear,  slip  and  partial 
shear  connection  significantly  affect  deflection  values 
and  should  therefore  be  considered  for  an  accurate 


evaluation  of  deflection. 


TABLE  6.1  SHEAR  FORM  FACTORS  (K  ,  )  FOR  COMPOSITE  BEAMS 
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TABLE  6.2  DEFLECTION  COMPARISON  OF  SIMPLY  SUPPORTED  COMPOSITE  BEAM  AT  SERVICE  LOAD 
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FIGURE  6.2  SHEAR  FORM  FACTOR  (K  )  VS  SLAB  WIDTH  (be)  FOR  COMPOSITE  BEAMS 
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FIGURE  6.5  LOAD  DEFLECTION  RELATIONSHIP  (SIMPLY  SUPPORTED  BEAM, 
FULL  AND  PARTIAL  SHEAR  CONNECTION,  CONCENTRATED  LOAD  AT  CENTER) 
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FIGURE  6.9  LOAD  DEFLECTION  RELATIONSHIP  (CONTINUOUS  BEAM, 
CONSTANT  'I',  FULL  AND  PARTIAL  SHEAR  CONNECTION) 
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FIGURE  6.11  LOAD  DEFLECTION  RELATIONSHIP  (CONTINUOUS  BEAM, 
VARIABLE  'I',  FULL  AND  PARTIAL  SHEAR  CONNECTION) 
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FIGURE  6.13  COMPARISON  OF  DEFLECTION  VALUES  (CONSTANT  AND  VARIABLE 
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FIGURE  6.14  COMPARISON  OF  SLIP  VALUES  (CONSTANT  AND  VARIABLE 
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FIGURE  6.15  COMPARISON  OF  SLIP  STRAIN  VALUES  (CONSTANT  AND  VARIABLE 
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figure  q  17  comparison  of  deflection  values  (shored  and  unshored  construction) 
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FIGURE  6.18  COMPARISON  OF  SLIP  VALUES  (SHORED  AND  UNSHORED  CONSTRUCTION) 
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FIGURE  6.19  COMPARISON  OF  SLIP  STRAIN  VALUES  (SHORED  AND  UNSHORED  CONSTRUCTION) 
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CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 


7 . 1  Summary 

This  study  has  examined  the  formulation  of 
governing  equations  for  composite  beams.  Solutions  of 
these  equations  have  emphasized  the  effects  of  such 
factors  as  partial  and  full  shear  connection,  shear 
deformation  and  linear  or  nonlinear  load-slip  character¬ 
istics  of  the  shear  connectors. 

1.  Classical  solution  techniques  were  examined  in  order 
to  provide  a  standard  comparison  against  which 
numerical  solutions  can  be  measured. 

2.  Programs  were  written  for  the  evaluation  of  classical 
solutions  for  a  simply  supported  beam  and  a  continuous 
beam  of  two  unequal  spans,  subjected  to  concentrated 
loads . 

3.  Results  obtained  from  numerical  technique  developed 
herein  were  compared  with  published  numerical  results 
of  Yam  and  Chapman  ( 31) ,  test  results  of  Chapman 

and  Balakrishnan  (  5  )  and  Hamada  and  Longworth  (  9  ) . 
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4.  The  solution  for  the  deflection  of  a  simply 
supported  composite  beam  with  uniform  load  was 
recast  in  terms  of  two  dimensionless  parameters, 

£  and  n •  This  permitted  a  plot  to  be  developed 

to  determine  the  effect  of  slip  which  is  applicable 
to  sections  of  any  simply  supported  composite  beam. 

5.  Shear  form  factors  were  computed  for  various  rolled 
steel  shapes  with  different  concrete  slab  widths 
and  thicknesses.  These  factors  can  be  used  to 
compute  shear  deflection. 

6.  The  effects  of  parameters  such  as  partial  and  full 
shear  connection,  shear  deformation,  linear  and 
nonlinear  load-slip  characteristics  on  deflection  and 
slip  for  a  simply  supported  beam  have  been  examined 
for  two  loading  cases,  viz.  uniformly  distributed 
load  and  concentrated  load. 

7.  A  numerical  solution  technique  which  permits  the 
solution  for  slip  and  deflection  of  continuous  beams 
under  arbitrary  loading  was  developed  and  programmed 
This  program  is  capable  of  treating  the  effects  of 
slab  cracking  in  negative  moment  regions  and  non- 
linearities  in  the  shear  connector  load-slip 
relationship.  The  numerical  solution  technique  was 
checked  against  closed  form  solutions  and  other 
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numerical  solutions  and  results  were  compared  with 
selected  experimental  results.  The  program  was 
used  to  prepare  a  simple  design  chart  to  determine 
deflections  of  two-span  uniformly  loaded  continuous 
beams . 

8.  The  solution  for  deflection  due  to  shear  and  slip 

of  a  two  span  continuous  composite  beam  with  uniform 
load  was  recast  in  terms  of  the  same  two  dimensionless 
parameters  £  and  n  .  Plots  were  developed  for  beams 
with  uniform  and  variable  moment  of  inertia. 

9.  The  effects  of  various  design  parameters  on  deflections 
were  examined  for  particular  example  problems. 

10.  The  effect  of  shored  and  unshored  construction  on 
the  behavior  of  three  span  continuous  beams  was 
examined . 

7 . 2  Conclusions 

On  the  basis  of  this  investigation,  the 

following  conclusions  are  drawn : 

1.  There  is  good  agreement  between  the  results  of  the 
numerical  solution  and  the  classical  solution  for 
deflection  and  slip  for  simply  supported  and 
continuous  composite  beams.  However,  there  is  a 
significant  difference  in  the  slip  strain  at  the 
concentrated  load  locations  in  the  case  of 


continuous  beams. 
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2.  The  results  for  deflection  and  maximum  slip  using 
the  technique  developed  herein  are  in  good  agreement 
with  Yam  and  Chapman's  results,  except  for  the 
distribution  of  slip  along  the  length  of  the  beam, 
which  may  be  due  to  the  error  in  the  assumed  data. 

3.  The  results  of  the  numerical  solution  are  in  good 
agreement  with  the  test  results  of  Chapman  and 
Balakrishnan  in  the  study  of  linear  load-slip 
behavior.  However,  the  numerical  solution  results 
in  higher  deflection  values  than  test  results  for 
nonlinear  load-slip  behavior.  On  the  basis  of  the 
end  slip  characteristics,  the  load-slip  behavior  is 
linear.  The  results  obtained  with  the  conditions  of 
full  shear  connection  and  non-linear  load  slip 
behavior  underestimates  the  deflection,  whereas  the 
condition  of  partial  shear  connection  overestimates 
the  deflections  obtained  by  Hamada  and  Longworth  (  9  ) 
at  maximum  service  load,  even  though  Hamada  and 
Longworth' s  beam  was  designed  for  full  shear 
connection.  The  discrepancy  in  the  deflection  may 

be  due  to  the  assumed  connector  stiffness. 

4.  Deflection  due  to  slip,  in  the  case  of  simply 
supported  beams  subjected  to  uniformly  distributed 
load,  may  be  as  high  as  45%  of  the  bending  deflection, 
depending  on  the  stiffness  of  the  connector  and 
geometry  of  the  beam.  Deflections  due  to  shear  and 
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slip  can  be  computed  for  practical  design  purposes 
using  the  table  of  shear  form  factors  as  well  as 
the  plot  developed  in  terms  of  dimensionless 
parameters  £  and  p .  Partial  shear  connection 
increases  the  slip  and  reduces  the  effective  stiff¬ 
ness  of  the  beam,  thus  increasing  the  deflection. 
Therefore  the  effect  of  shear,  slip,  and  partial 
shear  connection  on  deflections  should  not  be 
neglected  in  design. 

5.  Using  the  plot  in  Fig.  6.8,  the  total  deflection 
due  to  shear  and  slip  can  be  computed  for  a  two 
span  continuous  beam  with  uniform  load.  Total 
deflection  can  be  as  high  as  three  times  the  bending 
deflection,  depending  on  stiffness  of  shear  connector 
and  beam  geometry.  The  effects  of  the  various 
parameters  on  deflection  and  slip  of  continuous  beams 
are  similar  to  those  for  simply  supported  beams. 
However,  increases  in  deflection  and  slip  tend  to  be 
much  higher. 

6.  Concrete  cracking  in  negative  moment  regions  produces 
increased  deflection  and  slip  in  positive  moment 
regions.  The  increase  in  deflection  may  be  as  high 
as  15%  of  the  deflection  of  composite  beams  with 
uniform  moment  of  inertia.  Cracking  also  creates 
discontinuities  in  the  slip-strain  relationship. 
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7.  Unshored  construction  results  in  increase  in 
deflection,  decrease  in  slip  in  both  positive 
and  negative  moment  regions,  decrease  in  slip 
strain  as  well  as  in  introducing  discontinuities 
in  the  slip-strain  relationship  and  decrease  of 
slab  force  in  both  positive  and  negative  moment 
region . 
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CLOSED  FORM  SOLUTION  (TWO-SPAN  CONTINUOUS  BEAM) 


A.l  Classical  Solution 


unequal 
of  each 


where  u 
k 


Figure  A1  illustrat 
spans  and  with  equal 

span.  The  deflection 

,4  ,2 

d  u  ^  d  u 

4  K  2  ~  ~ 

dx^  dx^ 

=  deflection 
(ZEI +  EAz2) 


es  a  continuous  beam  with  two 


concentrated  loads  at  the  center 


equation  may  be  expressed  as 


w 


M 


IE  I 


EAZEI 


(Al) 


EA^EI 


P  =  K  =  shear  connector  modulus 
M  =  bending  moment 
w  =  uniform  load 
W  =  concentrated  load 


I  E I 

=  E1I1  + 

E2I2 

1 

E1A1  + 

E1A1 

EA  E1A1E2A2 


where  E^  and  E ^  are  moduli  of  elasticity  of  concrete  and 
steel,  respectively;  A^  and  A ^  are  areas  of  concrete  slab 
and  steel  section,  respectively;  and  1^  and  I ^  are  moments 
of  inertia  of  concrete  slab  and  steel  section,  respectively. 


The  geometry  of  the  beam  is  as  shown  in 
Fig.  Al.  The  equations  for  moment  for  the  four  segments 
of  the  beam  are  as  follows: 


M^X) 


o  <  x  < 


( A2 ) 
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WL. 

M2(X)  =  ~2 


W  ^  Mo\ 

2  +  L-JX 


L1 

T<x<L1 


(A3) 


WL, 

M3^>  =  -T 


w 

2 


Mo 


y 


L. 


o  <y  <- 


( A4 ) 


WL, 

M4(y)  =  V 


W  ^  Mo  \ 

2  +  L-Jy 


L2 

~T<y<h2 


(A5) 


where  Mo  is  the  redundant  moment  at  the  interior  support 
W  W 


FIGURE  Al.  TWO-SPAN  CONTINUOUS  COMPOSITE  BEAM 


Equation  Al  may  be  solved  in  each  of  the 
four  segments  of  the  beam  as  follows: 

U 

SEGMENT  1  o  <  x  < 

Homogeneous  Solution: 

u^  =  C2  +  C^x  +  sinh  /kx  +  cosh  /kx 


(A6) 
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Particular  Solution: 

x3  W  _  Mo 

P  6El  2  L1 

Complete  Solution: 

u(x)  =  +  C^x  +  sinh  /kx  +  cosh  /kx 

+  ZL__  W  _  Mo 

6ET  2  L1 

SEGMENT  2  L1 

- —  ~2  <  X  <  L1 

Homogeneous  Solution: 


u,  =  C-  +  C„  x  +  Cr,  sinh  /kx  +  co  cosh  /kx 

h  5  6  7  8 


Particular  Solution: 


u 

P 


WL1x 


4EI  6EI 


x3  l W  Mo 
2  -  Li, 


=  ( 

vi  \ 


Complete  Solution 


u(x)  =  Cr-  +  Cnx  +  sinh  /kx  +  CQ  cosh/kx 
00/  o 


WL1x 


+  _ ^ _ x/  /  W  +  Mo 

4ET  6ET  V2  L1 


Similarly  by  symmetry: 


L, 


SEGMENT  3 


o  <  y  < 


u(y)  =  Cg  +  C1Qy  +  C^  sinh  /ky  +  C19  cosh  /ky 


'll 


12 


( A7 ) 


( A8 ) 


(A9) 


(A10) 


(All) 


6EI 


(A12) 
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SEGMENT  4 


u(y)  =  C13  +  C14y  +  C15  sinh  /Ty  +  C  cosh  /ky 


+ 


WL2y 


_ _  y  j  w  +  Mo 

4ET  6ET  '  2  L2 


( A13 ) 


where  C. 

l 


i  -  1,  2,  . .  16  are  constants  of 

integration  for  a  particular  beam,  loading  condition  and 
geometry . 


For  any  segment  of  the  beam,  the  equations 
for  force,  slip,  and  slip  strain  are  given  by: 


Force : 

ZEI 

P  =  z 

d~u  M 

,  2  z 

dx 

(A14) 

Slip  : 

s  -  1 

y 

dP 

dx 

(A15) 

Slip  strain : 

,  ds 

dx 

( A16 ) 

The  slope  at  any  point  is  defined  as 


9  =  ^ 
dx 


(A17) 


In  order  to  solve  for  the  17  unknowns  (16  constants  of 
integration  and  1  redundant  moment),  it  is  necessary  to 
utilize  17  boundary  conditions:  six  related  to  deflection, 
five  related  to  force,  three  related  to  slip  and  three 
related  to  rotation.  These  boundary  conditions  are 
illustrated  in  FigureA2.  A  17  x  17  matrix  was  formulated 
using  the  above  conditions  applied  to  the  equations  for 
deflection,  force,  slip,  and  rotation.  The  matrix  was 
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coded  in  Fortran  and  solved  using  the  subroutine  LINV3F 
in  the  subroutine  library  IMSL. 

After  the  constants  of  integration  and 
redundant  moment  were  solved,  values  for  moment,  deflection, 
force,  slip,  and  slip  strain  were  computed  at  3-inch 
intervals  along  the  length  of  the  beam. 


W  W 


FIGURE  A2 .  BOUNDARY  CONDITIONS  FOR  CONTINUOUS 

COMPOSITE  BEAM 


A.  2 


Summary  of  Equations 


Deflection 


In  segment  1 


x  /  W  Mo 


-  C  +  C9x  +  C^sinh  /kx* +  C.cosh  /kx  +  . 

6EI  L1 


In  segment  2 

rj  _ _  WL1  x' 

u9  =  Cp.  +  C  x  +  C  sinh  / kx  +  C0cosh  /kx  +  - 

2  0  6  8  4EI 


( A18 ) 


-  ~(l +  it) 

6EI  I  1/ 


(A19) 


In  segment  3 


u. 


=  0^  +  C-,  ny  +  C11sinh  /ky  +  C.,^cosh  /Tcy 


1CK  '  —o'  ■  -16' 

y"  (  W  Mo" 


+ 


6EI 


2  L, 


( A20) 


In  segment  4 


U4  =  ^13  +  ^14y  +  ^15  s:i-n^  ^ ky  +  cosh  /  ky 


+ 


WL^ 

4ET 


W  +  Mo 
6EI  V  2  L2 


(A21) 


Slope 


In  segment  1 


A  . 


x  /W  Mo 


0  =  C  +  /kC  cosh  /kx  +  /kC  sinh  /kx  +  _ ,9  - 

6  2e I  o  -  L- 


(  A22 ) 
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In  segment  2 


WL.x 


-  Cg  +  /k*C7cosh  /kx  +  /"kC0sinh  /kx  +  — — 


'1 


8k 


2EI 


x /W  Mo 

2ET  V2  L1 


In  segment  3 


In  segment  4 


-  C-^o  +  /kC/  1  cosh  /ky  +/kC19sinh  /ky  + 


y 


li 


12‘ 


2EI 


04  =  +  /kC1  (.cosh  /ky  +  /kC1c;sinh  /ky 


15 


16 


+ 


wv 


y 


w  +  Mo 
2EI  2ET  V  2  L2/ 


Force 


In  segment  1 

P  =  ^-^{kCqsinh  /kx  + 

1  z  y  6  2EIV 


Mo\\_  x/W  _  Mo' 
Ljj-  z\ 2  -  L. 


In  segment  2 

P9  =  ^  E ^  /kC^sinh  /  kx  +  kCoCOsh  /lex  +  — -  - 
Z  Z  {  f  6  OFT 


x_/V 

eT 


1/WL1  /W  ,  Mo\ 
z(  2  (2  X 


In  segment  3 
EEI 


\ 

< 

/ 


P3  +  z 


{kCllsinh  /ky  +  /;(|  + 


y  f  w 

Z  V  2 


Mo 

Lr 


In  segment  4 
ZEI 


P4  =  z 


jkC^gSinh  /ky  +  kC^gCOsh  / ky  + 


WL. 


2EI 


y_/W  +  Mo 
ET\2  L2 


-  i{ 


WL2  (  W  ^  Mo^ 
2-  -(  2  + 


(  A23 ) 

W  Mo 
2  L2 
(  A24 ) 

(  A25 ) 


( A26 ) 


( A27) 


(  A28 ) 


( A29  ) 
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Slip 


In  segment  1 


1 


-  z  k  ^  Crj  cosh  /ky  +  k3//2C^  sinh  / kx 

+  W  _  Mo _ l/w  _  Mo\ 

El  2  Li  yz\ 2  Li  J 


In  segment  2 

s2  =  ^  {  k3/2Cy  cosh  /ky  +  k3/2C.  sinh  /kFF 


_1_  W  +  Mo  +  j/w  Mo 

IT  2  Li  yz\2  Li 


In  segment  3 

=  Jli  {  k3/2C„  cosh  /k?  + 
3  y  z  l  11  17 


1  W  Mo 


+  k3//2C^2  sinh 


El  2  L2 

l/W  Mo' 
z  V* 1 2  L2 


In  segment  4 

’4  =  /IT  {k3/2cosh  /ky  +  k3//2C16  sinh  /ky 


1/W  Mo  \  \  1/W  Mo 

=V2  L2H  zy2  L2 


Slip  Strain 

In  segment  1 

s'  =  i  k2C0  sinh  /kx  +  k2C/1  cosh  /kx  1 

1  y  z  l  3  4  ) 

In  segment  2 

s'  =  {  k 2 C rr  sinh  /kx  +  k2CQ  cosh  /kx  ( 

2  y  z  l  7  o  ’ 


In  segment  3 

Z  El  |  ,_2 


3  yz  ^11  sinh  /ky  +  k^C 


12 


( A30 ) 


(  A  3 1 ) 


( A32 ) 


( A33 ) 


( A34) 


( A35) 


cosh 


( A36) 
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In  segment  4 


s ' 


4 


ZEI  f 

y  z  l 


_  2 

sinh  / ky  + 


cosh 


( A37 ) 
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A . 3  Equations  Used  in  Solution  of  Constants  of 
Integration  and  Redundant  Moment 


(1)  c±  +  c4  =  0 


(2)  C9  +  C12  =  ° 


(3)  C-  + 


(5) 


'1/2  6 


(4)  Cn  + 


13 


y/kL1 

/KT, 

+  sinh  — 

2  C3  +  cosh  2  C 

/kL-. 

-  sinh  — g 

—  Crj  -  cosh 

to 

h 

n 

00 

i 

tr* 

i— 1 

3 

24 

El 

/kL9 

/kr2 

C^q  +  sinh 

2  C  + 

2  U11 

cosh  — 2 — 

/k'L2 

/kL' 

sinh  — ^ — - 

C,  _  -  cosh  — 
15 

2  C 

2  16 

WL1 

3 

C 


12 


6EI 


C,  +  L-Ce  +  sinh/kL1C7  +  cosh  /kL-.  CQ  -  — 
516  1/  -L  o  a 


Lh  2 Mo 


6EI 


WL. 


6EI 


(6) 


/kL/  /kL^  t/IcL^ 

ZEIk  sinh  ^  -  -  C„  +  cosh  — C,  -  sinh  -  Cr 


2  3 


2  4 


/kL, 


C8  = 


cosh 


2 


0 
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(7)  ZEIk  sinh 


^kL2  /kL2  Al2 

2  C11  +  cosh  2  C12  '  Sinh  “T"  C15 


/kll 


-  cosh 


2  C16  ~  0 


(8)  kC4  =  0 


(9)  kC12  =  0 


(10)  ZEIk  (sinh/kL1C7  +  cosh^in^Cg  -  sinh/kL9C 


15 


(11)  C2  +  /kcosh  —  2 


-  coshAL2C16)  =  0 
SkL1  /klT 


+  /ksinh 


C 

2  ^4 


/kL  /klT  WL  2 

-  /kcosh  — 9 —  C7  -  /kCQsinh  — » —  =  - — 

2  7  8  2  gEI 


(12) 


Ikl*  TkL 

CL  ~  +  /kCL  ..  cosh  — —  +  /kCL^sinh 


10 


'll “  2 

/kL, 


12' 


_  r 
2  14 


/kL  WL 

-  /kC1£-cosh  — ~ — -  -  /kc^^sinh  — ^ —  =  — — 
15  2  16  2 


8  El 


(13)  C~  +  /kC„cosh  /kL  +  v^c"0sinh  /kL-,  +  C1  . 

b  /  1  o  1  14 


Mo 


+  /kC^j-cosh  /~kL2  +  ZkC^gSinh  /kL2  -  ^  (L2  +  L1) 


2EI 


Z=  (L/  +  Lx2) 


4EI 


(14) 


3/2  /  /kL^ 

ZEIk  7  (  C3cosh 


/kL. 


2  4-  C4sinh  2" 


/kL. 

C7cosh  — ~ 


/kL. 

-  CgSinh  — 2~ 


1  -  ^  )  W 
El 
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(15) 


(16) 


(17) 


ZEIk3/2  (  C-L1cosh  — 


/kL, 

+  C12sinh 


-  C 


-  C-^sinh  — 
16  z 


=  W  1  - 


El 

EL 


3/2 


£EIk~'“  1  C^cosh  /kL^  +  Cgsinh  /kL^  +  C^5 


El 


+  C^sinh/kO  +  Mo  = —  +  t —  1  -  —  1  =  W 
lb  2]  L1  L2  \  El 


C13  +  Ci4L2  +  C15sinh/kL2  +  C16cosh  /kL2 

3 


-  Mo 


L2z 

6ET 


WL, 


6  El 


-  cosh 
15 


2 


cosh  /kL2 


Constants  Used  in  Closed  Form  Solution 


Loading : 

W  =  22.5  kips 

w  =  0.0  kips/in. 

Lengths : 

L^  =  180  in. 

L 2  =  120  in. 

Shear  Connector  Modulus:  y  =  K  =  150  kips/ 
Beam  Section  Properties: 


Therefore 

=  3000.0  ksi  (transformed) 

E2  =  30,000.0  ksi 

I  =  32.0  in.4 

I9  =  278.28  in.4 

=  24.0  in . 2 

A9  =  10.499  in.2 

z  =  8.12  in . 

9 

I El  =  9312870.0  kips-in.~ 

2 

El  =  23760159.0  kips-m. 

EA  =  219115.9  kips 

k  =  0.0017465 
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LINEAR  CLOSED  FORM  SOLUTION  FOR  A  SIMPLY 
SUPPORTED  BEAM  WITH  CONCENTRATED  LOAD  AT  CENTER 

B .  1  Classical  Solution 

Figure  B1  illustrates  beam  geometry,  loading 
pattern  and  boundary  conditions.  The  governing  equation 
for  slip  is  Equation  3.23, 

s"  -  a2s  =  -  BV  (Ba) 

The  homogeneous  solution  for  this  equation  is 


S^  =  A  sinhax  +  B  coshax 


(  B2  ) 


and  the  particular  solution  for  the  concentrated  load  at 


the  center 


s  V(X) 

a 


( B3 ) 


As  the  load  is  at  the  center  of  the  span 

V(x )  =  AP 


(B4) 


Therefore  the  general  equation  for  slip  can  be  expressed 


as 


s  =  A  sinhax  +  B  coshax  + 


3P 

o  2 

2a 


( B5  ) 


The  boundary  conditions  relating  to  slip  and  slip  strain 


are 


at  x  =  0 


M  =  0 

dx 


( B6  ) 


at  x  =  L/2 


s 


0 


(  B7 ) 
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Differentiating  Equation  B5  results  in  slip  strain  as 


^7  =  Aacoshax  +  Basinhax  (B8) 

By  substituting  boundary  condition  Equation  B6  into 
Equation  B8, 

A  =  0  ( B9 ) 


Substituting  boundary  condition  Equation  B7 
B5  gives 

st  /o  =  0  =  0  +  B  coshaL/2  +  - 

W  ^  n  Z 

2a 


Therefore 


D  8  P  ,  L 

B  =  -  — £  sech  a  2 

2  a 


By  substituting  constants  into  Equation  B5 


8  P  ^  L  ,  ,  8P 

s  =  -  2  sech  a  2  coshax  + 


2  a 


2a 


or 


s  =  ^  ^  (1  -  secha  ^  coshax) 
2  a 


and 


s '  = 


8  P 


L 


a 


2  ( -  a  secha  ^  sinhax) 


into  Equation 


(BIO) 


(Bll) 


( B12 ) 

(  B13 ) 

( B14 ) 


The  deflection  may  now  be  determined  by  integration  of 
Equation  B14  and  using  Eqaution  3.19(a).  Combining  these 
two  equations  yields 


_v”  = 


AYs’ 
c  c 

h 


+ 


M 


El 


(  B15  ) 


AY  t-j  q 

c  c  P  8 


2  oi 


,  .  L  .  ,  .  Px 

„  (—secha  77  sinhax)  +  tvft- 
z  z  zti  1 


t 


-v 


( B16 ) 
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Integrating  Equation  B16  yields 


P  BA  Y 

c  c  2 

-V  =  -  - ft  a  sechaL/2  coshax  Px  ^ 

2 1 1  az  a  4EI  L1 


(B17) 


P8A  Y  ,  T  /o  -v.  3 

t t  T7i  c  c  sechaL/2  sinhax  Px  r  „  ,  r  /P1on 

v  „T  2  5  +  T2ETT  +  Clx  +  C2  (B18> 

Zit  t 

The  boundary  conditions  relating  the  deflection  and  slope 


are 


at  x  =  0 


at  x  =  L/2 


v  =  0 


v'  =  0 


( B19 ) 
( B20 ) 


By  substituting  Equation  B20  into  Equation  B17  yields 

2 


P  BA  Y 

0  =  -  - -  ft  secha^  cosha-^  +  +  CL 

2Ita  t 


(B21) 


Therefore 


ci  = 


P8A  Y 
c  c 


PL' 


2 1  ta 


2  16E I 


( B22 ) 


Substituting  Equation  B19.  into  Equation  B18  yields 


C2  =  0 


(B23) 


Substituting  Equation  B18  yields 


P®AcYc  L  Px3  P6AcYc 

-  v  =  -  - o  (sechayr  sinhax)  +  ? 

2T  a  t  21^  a' 


t 


t 


PL 


16E I 


x 


(  B24 ) 


t 


Equations  B13,  B14  and  B24  give  the  solution  for  slip,  slip 
strain  and  deflection  of  a  simply  supported  beam  subjected 
to  a  concentrated  load  at  midspan.  At  x  =  L/2  Equation  B24 


. 
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yields  the  maximum  deflection  in  the  span  as 


or 


VL/2 

VL/2 


PL* 


48E I 


+ 


PL' 


48EI 


t 


P3A  Y 

CO,  i  p.  I— i 

- o  tanha—  - 

21  &  Z 
PBA  Y 

+  - — ^(tanha-^ 

2  ItaJ 


PBA  Y 
c  c 

4It“2 


Lx 

a2} 


( B25 ) 


( B26 ) 


Recognizing  that  the  first  term  is  the 
beam  with  complete  interaction  Equation 

aS  6  „  24BEA  Y 

f  c  c  ,  L  aLN 

—  =  1  +  3  3  (tanha5  -  -g-) 

o  La 


deflection  of  a 
B26  may  be  written 

( B27 ) 


Recalling  the  definition  of  a  and  B  from  Equations  3.21 
and  3.22,  Equation  B27  may  be  written  in  terms  of  two 
nondimensional  ratios  as 


6 - 1  + 

o 


(?n) 


245  ,  .  /?n 

—3/2  (tanh  7T 


or 


o 


=  1  +  24 _  ftanh  4a 

1  .1/2  3/2  (-ta  2 

£  '  n  ' 


The  definition  of  £  and  q  are  given 


(  B28 ) 

^ )  ( B29 ) 

in  Equations  6.4  and 


6.5,  respectively. 
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B . 2  Evaluation  of  Shear  Deflections  for  a  Simply 

Supported  Composite  Beam  Subjected  to  a  Concentrated 

Load  at  the  Center 


The  shear  strain  in  a  beam  may  be 


expressed  as  K  V 

Y  _  ,  _  sh 

v  s(x)  GA^ 


in  which  v  is  the  shear  deflection  and 
s 


( B30 ) 


( B31 ) 


is  a  property  of  the  cross  section.  For  a  simply  supported 
beam  subjected  to  concentrated  load  at  midspan, 


6 

s 


Shear  deflection  can  therefore 
is  known. 


be  evaluated 


Vdx  ( B32 ) 

provided  K  , 
^  sh 


The  integral  on  the  right  hand  side  of  the 
equation  is  the  moment  at  midspan.  The  ratio  of  shear 
deflection  to  bending  deflection  can  be  expressed  as 


K 


sh  PL 


48EI 


t 


or 


o 

6s 

6 

o 


GAt  4 


PL' 


=  12  K 


E 


t 


sh  G  AtL2 


By  adding  Equations  B29  and  B34, 


+ 


6  +  6  +  6 
o  si  s 


o 


o 


o 


( B33 ) 

( B34 ) 


( B35 ) 
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Equation  B35  is  the  ratio  of  total  deflection  to  bending 
deflection  of  a  simply  supported  composite  beam  subjected 
to  concentrated  load  at  midspan. 
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EVALUATION  OF  SERVICE  LOADS  OF 
COMPOSITE  BEAMS  FOR  TEST  PROBLEMS 

C . 1  Introduction 

Most  experimental  research  in  the  behavior  of 
composite  beams  is  concerned  with  ultimate  loading  con¬ 
ditions.  However,  this  present  study  is  concerned  with 
behavior  under  service  load  conditions.  The  following 
procedure  is  followed  to  evaluate  service  loads  for 
experimental  beams  so  that  results  from  this  analysis 
can  be  compared  with  test  results  reported  by  other 
researchers . 

Under  service  load  conditions  the  behavior  of 
composite  beams  is  governed  by  the  geometry  of  the 
composite  section  and  stiffness.  Under  ultimate  load, 
however,  the  behavior  is  a  function  of  geometry  of 
composite  section  and  material  strength,  that  is,  yield 
strength  of  steel  and  ultimate  strength  of  concrete  and 
shear  connectors. 

In  Chapter  V  the  experimental  results  obtained 
by  Chapman  and  Balakrishnan ,  and  Hamada  and  Longworth 
are  compared  with  analytical  results  obtained  from  this 
study  under  service  load.  The  properties  of  these 
beams  used  in  the  numerical  analysis  are  computed  in 


this  Appendix. 


C.  2 


Simply  Supported  Beam  with  Concentrated  Load 

at  Center  (Chapman  &  Balakr ishnan ) 


Figure  Cl  illustrates  the  beam  and  its  cross  section. 

Beam  Geometry: 

Length  =  L  =  18.0  ft. 

Width  of  slab  =  be  =  48.0  in. 

Thickness  of  slab  =  tc  =  6.00  in. 

Steel  section  BSB12  x  6  x  44# 

Depth  of  steel  beam  =  ds  =  12.0  in. 

Area  of  steel  beam  =  A  =  13.0  sq .  in. 

s 

Area  of  web  =  A  =  4226  sq .  in. 

w  ^ 

Thickness  of  flange  =  tf  =  0.717  in. 

Shear  Connectors: 

Size  of  shear  connector  3/4”  x  4” 

Spacing  of  shear  connector  =  14.88  in. 

Total  number  of  shear  connectors  =  N  =  32 

Material  Properties: 

Yield  strength  of  flange  =  f  =  34.832  ksi 

Yield  strength  of  web  =  f  =  38.8192  ksi 

yw 

3 

Elastic  modulus  of  steel  =  E  =  31.6  x  10  ksi 

c 

Concrete  strength  (28  day  cylinder)  =  f^  =  3440  ps 
Ultimate  load  capacity  of  connector  =  qr=  28  kips 
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Ultimate  capacity  of  concrete  slab  =  C  =  0.85  f^  bctc 

=  0.85  x  3.44  x  48.0  x  6.0  =  842.11  kips 

Total  capacity  of  shear  connector  =  =  Nqr/2 

=  28.0  x  16.0  =  448.0  kips 

Ultimate  capacity  of  steel  section  =  T  =  A  f 

s  y 

=  A  f  x  f  f  +  A  w  x  f  w  =  469.67  kips 
s  y  s  y  H 


Therefore  T  >  Q 


u 


T  <  C 


C>QU 

As  the  connector  capacity  is  less  than 

ultimate  force  that  can  be  developed  by  the  steel 

section  or  the  concrete  slab,  and  concrete  slab  strength 

is  greater  than  the  force  that  can  be  developed  by  steel 

section,  the  governing  criteria  in  determining  the 

location  of  the  neutral  axis  at  load  is  Q  .  Therefore 

u 

it  can  be  concluded  that  there  is  partial  shear  connection 
in  this  case.  To  achieve  full  shear  connection  the 
connectors  must  be  designed  to  transfer  the  ultimate 
force  that  can  be  developed  by  the  steel  section.  The 
degree  of  shear  connection  for  this  beam  is 

Q, 


Partial  shear  connection  (P.S.C.)  - 


u 

T 


448.00 

469.67 


x  100  =  95.38 


Due  to  this  partial  shear  connection,  the 
force  Qu  can  be  developed  in  the  steel  section  and 
tensile  force  must  be  equal  to  compressive  force. 


Force  that  can  be  developed  by  bottom  flange 
F1  =  Af  x  f  f  =  4.34832  x  34.832  =  151.46  kips 

Force  that  can  be  developed  by  web 

F2  =  Aw  x  f  =  4.2264  x  38.8192  =  164.06  kips 

The  amount  of  force  that  should  be  developed  in  top 
flange  to  utilize  the  capacity  of  shear  connection 

F3  =  Qu  -  Fp  -  F2  =  448.0  -  151.46  -  164.06  =  132.48  ki 

Depth  of  top  flange  below  ultimate  strength  neutral  axis 
dp  =  132.48/6.1183  x  34.832  =  0.62164  in. 

Depth  of  top  flange  above  ultimate  strength  neutral 
axis  under  compression 

d2  =  0.717  -  0.62164  =  0.9536  in. 

Compressive  force  that  can  be  developed  in  this  portion 
of  top  flange 

f4  =  0.09536  x  6.1183  x  34.832  =  20.322  kips 

The  magnitude  of  compressive  force  that  must  be  developed 

in  the  concrete  slab  for  equilibrium 

F  =  q  -  F .  =  448.0  -  20.322  =  427.678  kips 
5  u  4 

The  depth  of  stress  block  in  the  concrete  slab  under 
compression 

a  =  427.678/0.85  f'bc  =  3.0472  in. 

(3 
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Therefore  location  of  ultimate  strength  neutral  axis 
CG  =  (tf  +  dw  +  d1)  =  11.905  in. 


The  moment  (neglecting  strain  hardening)  that  can  be 
developed  in  the  composite  section  before  the  failure 
of  the  shear  connection 


Mp  =  F1(CG  -  V2)  +  F2(CG  '  tf  -  V2)  +  F3dl/2 

+  F4d2/3  +  F5(t  +  d2  -  a/2) 

=  4715  in.  kips  =  393  ft.  kips 

Mp  =  Moment  due  to  self  wt  +  Moment  due  to  applied  load 
Self  wt  =  44  +  -6~X-  i44X  150  =  344  lbs. 


Load  factor  for  dead  load  =  1.25 


Load  factor  for  applied  load  =1.5 


1.25  x  0.344  x  (118)  ^  1 . 5P  x  18 

M  =  - - - - - —  + 


P 


8 


P  =  (393.0  -  17.415)  x^=  83.463  kips 

U  lo 


( service ) 


83.463 

1.5 


=  55.65  kips 


THe  ultimate  load  given  by  Chapman  and  Balakrishnan 

P  =  96.95  kips  and  the  load  factor  =  1.75 
u  1 


p 

( service ) 


96.95 

1.75 


55.4  kips 


Therefore  it  is  evident  that  service  load  obtained  by 
the  above  procedure  is  accurate  to  the  degree  required. 
Figures  El  to  E3  illustrate  the  beam  geometry  and  the 


loading  condition. 
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C.3  Two-Span  Continuous  Composite  Beam  Subjected  to 
Equal  Concentrated  Loads  at  the  Center  of  Each 
Span  (Hamada  &  Longworth) 

Figures  C2  and  C3  illustrate  the  section  and 
beam  geometry  of  the  test  specimen  and  Figure  C4 
illustrates  the  idealized  stress  condition  at  ultimate 
moment  in  positive  moment  region.  Figures  C5  and  C8 
illustrate  the  stress  condition  in  the  negative  moment 
region . 

Beam  Geometry: 

Length  =  L  =  12.0  ft. 

Width  of  slab  =  be  =  48.0  in. 

Thickness  of  slab  =  tc  =  4.00  in. 

Steel  section  =  W12  x  31 

Depth  of  steel  section  =  ds  =  12.09  in. 

Depth  of  web  =  d  =  11.16  in. 

^  w 

Width  of  web  =  w  =  0.265  in. 

w 

Area  of  web  =  A  =  2.9574  sq .  in. 

w 

Area  of  flange  =  A^  =  3.0863  sq.  in. 

Thickness  of  flange  =  tf  =  0.465  in. 

Width  of  flange  =  b^  =  6.6372  in. 

Total  area  of  steel  section  =  A  =  9.13  sq .  in. 

o 

Area  of  reinforcement  in  negative  moment  region 

=  A  =  1 . 6  sq .  in . 
s  r 
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Shear  Connectors: 

Size  of  connector  3/4"  x  3" 

Number  of  connectors  between  zero  and  maximum 

moment  =  N  =  16 

Capacity  of  a  connector  =  qr  =  28.716  kips 


Material  Property: 

Yield  strength  of 
Yield  strength  of 
Elastic  modulus  of 
Ultimate  strength 
Elastic  modulus  of 
Yield  strength  of 


flange  =  F 

yt 

web  =  F  „  = 

yf 

steel  =  Es 
of  concrete 
concrete  = 
reinforcing 


=  40.5  ksi 
46.9  ksi 

=  30.550  x  103  ksi 

=  f'  =  5577  psi 
c 

E  =  4.527  x  103  ksi 
c 

steel  =  F  -  50.3  ksi 
ys 


In  positive  moment  region 

Ultimate  strength  of  slab 

C  =  0.85  f '  bet  c  = 
c 


The  ultimate  strength  of 


T 


A  f 


s  y 


is 

910 . 1664  kips 

steel  section  is 

+  A  f  =  387.51 
sw  yw 


kips 


The  ultimate  capacity  of  shear  connectors 

Q  =  N.qr  =  16  x  28.716  =  459.456  kips 
u 

Therefore  Q  >  T 
u 

T  <  C 
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As  T  and  C  T  full  shear  connection  can 

be  achieved.  As  the  concrete  slab  ultimate  strength  is 
greater  than  ultimate  capacity  of  the  steel  section,  the 
ultimate  strength  neutral  axis  will  be  in  the  concrete 
slab.  The  depth  of  the  stress  block  in  the  concrete 
slab  under  compression  is 


T 

a  0.85  fr  be 
c 

Therefore 


2 . 0192  in . 


CG  -  ds  +  t  -  a  =  14.0708  in.  from  the  bottom  flange 


The  ultimate  moment  that  can  be  developed 
by  the  composite  section  in  the  positive  moment  region 
is 

Mul  =  T  (CG  -  ds/2  +  a/2)  =  3501.3078  in.  kips 

=  291.775  ft.  kips 

In  negative  moment  region: 

In  the  negative  moment  region  as  illustrated 

in  Figure  C4 ,  the  concrete  slab  is  not  effective  as  it 

is  subjected  to  tension,  therefore  the  moment  has  to  be 

resisted  by  the  steel  section  and  the  reinforcement.  From 

Figure  C6  the  equilibrium  equation  can  be  expressed  as 

CS  +  T  =  TS 
r 

but  CS  +  TS  =  A  f 

s  y 

T  =  1.6  x  50.3  =  80.48  kips 
r 


CS  +  TS  =  387.51  kips 
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Therefore  CS  =  (387.51  -  80.48)/2  =  153.26  kips 
TS  =  387.51  -  153.26  =  234.25  kips 

The  bottom  flange  can  develop 

F1  =  Af  x  f  f  =  3.0863  x  40.5  =  125  kips 

The  web  below  the  ultimate  strength  neutral  axis  must 
resist 

F2  =  234.25  -  125.0  -  109.25  kips 

The  depth  of  web  below  ultimate  strength  neutral  axis  is 

d  ,  =  109.25/46.9  x  .265  =  8.7902  in. 
wl 

Therefore 

CG  =  8.7902  +  .465  =  9.2552  in. 

The  depth  of  web  above  the  ultimate  strength  neutral 
axis  is 

d  0  =  (11.16  -  8.7902  =  2.3698  in. 
w2 

Area  of  web  above  neutral  axis  =  2.3698  x  .265 

=  0 . 628  sq .  i . n . 

The  force  that  can  develop  in  web  is 

F3  =  0.628  x  46.9  =  29.453  kips 

Force  that  will  develop  in  top  flange 

F4  =  3.0863  x  40.5  =  125  kips 

The  center  of  reinforcement  is  3  in.  from  the  bottom 
of  the  slab,  and  the  force  in  the  reinforcement  is 
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Fj_  =  80 . 48  kips 


Therefore  the  moment  that  can  be  resisted  by  the  section 
in  the  negative  moment  region  is 

M  0  =  F1  (CG  -  tf/2)  +  Fn  (CG  -  tf  -  d  1/2) 

+  F3  (Dw2/2)  +  F4  <dw2  +  tf/2) 

+  F,_  (ds  +  tc  -  CG  -  1.0) 

=  2437.773  in.  kips  =  203.147  ft.  kips 

^i]  9  2  03  147 

Total  M  =  M  ,  + — ^  =  291.766  +  -  l-1  *  •  =  393.35  ft.  kips 

p  ul  2  2 

WL^ 

Moment  due  to  self  wt  =  -  Q--  =  5.1975  ft.  kips 


M 

P 


M(due  self  wt )  +  M 


WL 

8 


+ 


u 


(due  to  applied  load) 
L 


P  L 

-^-=  393.35  -  5.1975  =  388.1525  ft.  kips 
4 

P  =  129 . 384  kips 
u 

Load  factor  =  1.5 

129  384 

Service  load  P  (service)  =  — -  =  86.256  kips 

_L  .  O 
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FIGURE  Cl.  SIMPLY  SUPPORTED  COMPOSITE  BEAM 
SUBJECTED  TO  CONCENTRATED  LOAD  AT  MIDSPAN 
(CHAPMAN  &  BALAKRISHNAN) 
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FIGURE  C2.  DETAILS  OF  TEST  SPECIMEN 

(HAMADA  &  LONGWORTH) 


181 


CQ 

< 

2 

CO 

2 

Fh 

2 

2 

o 

2 

o 

^  2 
O 


2 

Eh 

2 

O 

o 

f- — I 

o 

2 

Q 

2 

< 

< 

Q 

< 

g 

C 

2 


< 

2 

2 

Eh 

CO 

2 

2 

a 

2 

o 

2 

2 

Eh 

2 

O 

t— i 

U 

2 

CO 


♦ 

2 

O 

2 

2 

D 

O 

t— ( 
2 


.  85f '  0 . 85  f 


182 


2 

2 

S 

o 

s 

w 

c-1 

< 


J 

P 

c-1  O 
<  2 


o 


w 

2 

d 


CO  G 

2:  2; 
O  W 
h-i  G 

G  pq 
Q  > 

2  t-H 

O  Eh 

O  H-t 

CO 
CO  O 
CO  G 
W 

2  2 
Eh  HH 

CO 

Q 

W 

CS3 

HH 

G 

< 

H 

Q 


Eh 

2 

W 

s 

o 


o 

w 

2 

P 

o 

I— 1 
2 


183 


2 

o 

a 

t-H 

2 

a 

i—i 

p 

p 

03 

p 

H 

P 

HH 

2: 

>H 

p 

03 

O 

W 

s 

P 

P 

p 

< 

> 

(-H 

P 

P 

< 

s_x 

a 

p 

2; 

< 

2 

a 

i— i 

2 

i — i 

2 

Q 

O 

P 

l-H 

P 

P 

h- 1 

P 

P 

a 

HH 

p 

03 

aj 

P 

O 

CO 

P 

l-H 

p 

Q 

CQ 

CO 

CO 

aj 

p 

03 

P 

CO 

LO 

u 

p 

03 

P 

a 

t— t 

p 

184 


'XD 


>> 

4-1 


2: 

w 

s 

o 


w 

H 

< 


H 


H 

<C 

CO 

o  o 
>-<  ^ 
E-«  t-t 
D  Q 
CQ  Co 

HH  « 

C£  ca 

H 

co  w 

>— i  > 

Q  hh 

E~< 

CO  < 
CO  o 
H  W 
ca  S 
E-< 

CO 

Q 

W 

CO 

I— I 

J 

< 

w 

Q 


CD 

U 

w 

C£ 

CD 

c 

KH 

ptt 


CO 

3 


APPENDIX  D 


186 


EVALUATION  OF  SHEAR  CONNECTOR  STIFFNESS 
FOR  VARIOUS  TYPES  OF  SHEAR  CONNECTORS 


D.l  Evaluation  of  Shear  Connector  Stiffness 


In  Chapter  III  the  shear  connector  stiffness 
K  used  in  analysis  has  been  defined  as  load  per  connector 
per  unit  slip.  In  the  numerical  solution  shear  connection 
stiffness  has  been  divided  by  the  connector  spacing  as 
the  computer  program  is  based  on  the  numerical  technique. 

For  linear  load  slip  behavior 

Q  =  Ks  (Dl) 


fJQ  =  K 

ds 


For  nonlinear  load  slip  behavior 


(D2) 


Q  =  a(l  -  e  bs)  ( D3 ) 

§  =  ab  =  K  ( D4 ) 

Figures  Dl  to  D5  illustrate  load  slip  character¬ 
istics  for  various  types  of  shear  connectors  (Chapman  and 
Balakr ishnan ) .  For  linear  behavior  the  shear  connector 
stiffness  has  been  evaluated  as  the  slope  of  the  load- 
slip  curve.  For  nonlinear  behavior  the  value  of  a  and  b 
has  been  evaluated  using  the  procedure  as  described  in 
Section  2.1.3  (Yam  and  Chapman). 

In  most  of  the  examples  used  for  verification 
of  the  numerical  analysis,  difficulty  was  encountered  in 
obtaining  a  shear  connector  stiffness.  The  example 
problems,  which  are  compared  with  the  closed  form 
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solution,  are  designed  for  full  shear  connection  and 
3/4"  x  4"  stud  shear  connectors  are  used  to  achieve 
this.  The  load-slip  relationship  for  this  connector 
is  available  from  push-out  test  results  (Chapman  & 
Balakrishnan ) . 

In  the  example  based  on  Chapman  &  Balakrishnan ' s 
beam,  the  shear  connector  stiffness  was  obtained  from 
their  load-slip  data.  As  the  ultimate  strength  of 

concrete  in  the  push-out  test  and  the  test  beam  are 
different,  a  linear  interpolation  was  used  to  obtain  the 
shear  connector  stiffness.  The  difference  in  results 
may  be  attributed  to  this  interpolation. 

In  the  case  of  Hamada  and  Longworth’s  beam, 

CBI ,  no  push-out  test  results  were  available.  To  obtain 
a  shear  connector  stiffness  for  this  example,  several 
examples  were  selected  with  various  shear  connector 
stiffnesses.  After  analysing  these  examples,  the  end 
slip  from  analysis  was  plotted  against  shear  connector 
stiffness.  With  the  end  slip  for  Hamada  and  Long- 
worth's  beam,  a  shear  connector  stiffness  was  interpolated 
from  this  plot.  Figure  D6  illustrates  the  relationship 
between  end  slip  and  shear  connector  stiffness. 
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Push-out  test  results  for  various  stud  shear 
connectors  are  illustrated  in  Fig.  D1  to  D5.  The 
ultimate  load  for  these  connectors  are  also  obtained 
from  test  results.  Two  slip  values  are  selected  as 
=  0.01  and  y^  =  0.02.  Expressing  them  in  terms  of 
ratio  of  Q/Qu  (  ) 


a  = 


Qp  2 

^U 


2Q. 


Q 


u 


Q. 

Q* 


u 


(D5) 


b  = 


Q- 

Q 


Y 


£n 


u 


1 


q2  qx 


Q„  Q 


u 


u 


(D6) 


For  3/4"  x  4"  stud  connector 


Q  =  12.5  tons 


For  y  ^ 


Q- 

Q 


=  0.35555 


u 


For  Y  Q9 

z  =  0.51666 


Therefore  a  =  8.1j26  tons/in.  or  18.20  kips/in. 
b  =  79.1578 
K  =  1440.67  kips/in. 


For  other  stud  shear  connectors,  the  values  of  a,  b, 
and  K  are  computed  as  per  the  above  procedure  and  are 
tabulated  in  Table  D1 . 


In  the  analysis  the  shear  connector  stiffness 
was  divided  by  the  connector  spacing  in  a  particular 
beam  to  obtain  the  shear  connector  stiffness  per  unit 


length  of  beam. 
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FIGURE  Dl.  LOAD-SLIP  RELATIONSHIP  FOR  PUSH-OUT  SPECIMEN 
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FIGURE  D2 .  LOAD-SLIP  RELATIONSHIP  FOR  PUSH-OUT  SPECIMEN 
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APPENDIX  E 
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SHEAR  STRESS  DISTRIBUTION  IN  COMPOSITE  BEAM  SECTIONS 

E . 1  Evaluation  of  Shear  Stress 

In  a  composite  beam  the  neutral  axis  may  be 
either  in  the  concrete  slab,  in  the  top  flange  of  steel 
section,  or  in  the  web  of  steel  section.  For  these  three 
locations,  the  variation  in  statical  moment  of  area  and  the 
shear  stress  distribution  are  illustrated  in  Figs.  El  to  E3. 
The  shear  deflection  factor  for  three  common  steel  sections 
used  in  composite  beams  for  different  slab  widths  and  slab 
thicknesses  are  shown  in  Fig.  E4  to  E6 . 

The  shear  stress  distribution  is  based  on  the 
following  assumptions: 

(1)  There  is  linear  stress-strain  behavior  under 
service  load  conditions. 

(2)  Concrete  in  tension  is  neglected. 

It  can  be  observed  from  Fig.  E2  that  the 
maximum  shear  stress  occurs  at  the  neutral  axis  when  the 
neutral  axis  is  in  the  web,  whereas  in  the  other  two  cases 
(Fig.  El  and  Fig.  E3),  where  the  neutral  axis  is  in  the 
slab  or  in  the  top  flange,  the  maximum  shear  stress  occurs 
at  the  junction  of  top  flange  and  web.  Therefore  it  may  be 
concluded  that  the  major  portion  of  the  shear  force  in 
composite  beam  is  resisted  by  the  web  of  steel  section. 
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In  case  of  composite  beams  due  to  the 
presence  of  concrete  slab  on  the  top  of  steel  section, 
the  shear  stress  at  the  top  of  the  top  flange  of  steel 
section  is  now  zero,  therefore  the  shear  form  factors 
for  composite  beams  are  much  higher  than  that  of  wide 
flange  beams. 

Shear  form  factors  for  composite  beams  with 
steel  sections  (W16  x  50,  W16  x  26,  W14  x  53,  W14  x  22, 

W12  x  190,  W12  x  14)  for  various  slab  widths  are 
illustrated  in  Fig.  E4,  E5,  E6 ,  respectively,  for  concrete 
slab  thickness  6"  and  4". 
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FIGURE  El.  STATICAL  MOMENT  AND  SHEAR  STRESS 
DISTRIBUTION  IN  COMPOSITE  SECTION 
(NEUTRAL  AXIS  IN  SLAB) 


200 


SHEAR  CONNECTORS 


co 


o> 

o 

CM 


XX 


CONCRETE  SLAB 


NEUTRAL 


STEEL  SECTION 
( W12  x  31) 


AXIS 


(a)  SECTION  OF  COMPOSITE  BEAM 


(b)  STATICAL  MOMENT 
VARIATION 


(c)  SHEAR  STRESS 
DISTRIBUTION 


FIGURE  E2 .  STATICAL  MOMENT  AND  SHEAR  STRESS 
DISTRIBUTION  IN  COMPOSITE  SECTION 
(NEUTRAL  AXIS  IN  WEB) 
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COMPUTER  PROGRAM 


F.  1 


Introduction 


The  following  assumptions  are  made  in  the 
computation  of  slip,  shear  flow,  force,  slip  strain  and 
deflection ; 

1.  The  composite  beam  is  divided  into  a  number 
of  elements  of  constant  length,  thus  estab¬ 
lishing  a  number  of  equally  spaced  nodal 
points . 

2.  Discontinuity  of  strain  is  assumed  at  the 
interface  of  the  elements  due  to  the  presence 
of  slip  strain. 


O 

O 


Concrete  resists  no  tension  in  the  negative 
moment  regions. 

Reinforcement  exists  at  only  one  level  in 
the  slab. 

In  the  analysis  which  accounts  for  different 
section  properties  in  positive  and  negative 
moment  regions,  properties  are  constant 
throughout  each  region  and  equal  to  those 
associated  with  the  transformed  section 
(determined  for  full  interaction)  in  the 
respective  regions. 
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6.  Due  to  the  abrupt  change  in  section  properties 
at  the  junction  between  the  positive  and 
negative  regions,  the  slip  strain  is  dis¬ 
continuous  at  the  inflection  points. 

The  following  steps  are  required  in  the 
computations  for  slip,  slip  strain  and  deformation  of  a 
continuous  beam  for  a  given  load: 

1.  Initialization 

a.  Determine  the  neutral  axis  and  section 
properties . 

b.  Compute  the  stress  resultants  for  unit 
reactions  at  the  interior  supports. 

c.  Compute  the  stress  resultant  for  the 
external  load. 

2.  Determine  the  flexibility  influence  coefficients 
and  lack  of  compatibility  displacements. 

For  each  of  the  unit  reactions  in  1(b),  and 
for  each  of  the  stress  resultants  obtained  by 
combining  the  reactions  and  the  external  loads: 

a.  Compute  the  slip  at  each  nodal  point 
using  the  Runge-Kutta  method  of  numerical 
integration . 

b.  Iterate  the  slip  computation  until  the 
correct  slip  has  been  determined. 
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c.  Compute  slip  strain,  curvature  and 
deflection . 

d.  Enter  the  appropriate  deflections  in 
the  equations  of  consistent  deformation. 

3.  Solve  for  corrections  to  reactions. 

a.  Determine  corrections  to  redundant 
reactions  at  the  interior  supports. 

b.  Determine  the  location  of  the  inflection 
points  and  alter  the  length  of  beam  with 
negative  section  properties. 

c.  Check  the  deflection  at  the  supports 
computed  in  step  2(d). 

The  computer  program  was  written  in  Fortran 
IV  and  computations  were  carried  out  on  the  AMDHAL  470/v 
at  the  University  of  Alberta  computing  services.  The  flow 
chart  in  Fig.  B1  outlines  the  sequence  of  computation  required 
for  the  analysis  of  a  continuous  beam.  The  program  consists 
of  eight  subroutines  and  one  function. 

F . 2  Description  of  Program 

MAIN  PROGRAM:  The  main  program  reads  and  writes  the  beam 
dimensions,  material  properties  and  loading  con¬ 
ditions.  It  calls  the  following  subroutines  in 
sequence  to  compute  stress  resultants,  slip  and 


deflection . 
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PROP .  : 


REGULA 


CONC.  : 


BMSU : 


MASTER 


This  subroutine  computes  the  position  of  the  neutral 
axis  and  the  section  properties.  When  necessary  this 
subroutine  calls  subroutine  "REGULA"  to  compute  an 
average  width  so  that  average  section  properties  are 
computed  on  this  basis. 

This  subroutine  computes  average  width  of  concrete 
slab  by  using  the  section  properties  of  the  positive 
and  negative  moment  regions  by  using  the  REGULA- 
FALSI  method. 

This  subroutine  computes  stress  resultants  due  to 
all  concentrated  loads.  The  shears  computed  for  each 
element  are  stored  in  two  vectors,  identified  as  the 
left  shear  and  right  shear  of  the  element,  and  bending 
moment  at  the  nodal  point. 

Computes  stress  resultants  due  to  uniform  load  and 
adds  the  stress  resultants  to  those  computed  by 
"CONC"  for  external  concentrated  loads. 

This  subroutine  calls  two  subroutines  and  one 
function.  First  it  calls  "RUNGA"  to  compute  slip, 
shear  flow  and  axial  force.  A  correction  for  slip 
is  used  if  the  axial  force  computed  at  the  last  nodal 
point  is  not  zero.  For  variation  in  section  proper¬ 
ties,  this  subroutine  selects  the  section  properties 
depending  upon  the  sign  of  the  bending  moment  at  the 
nodal  point.  The  subroutine  also  computes  slip 
strain  and  curvature.  At  the  inflection  points  the 
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discontinuous  slip  strains  are  computed  at  the 
left  and  right  side  of  the  point.  The  subroutine 
calls  "DEFLEX"  to  compute  deflection  at  each  nodal 
point  assuming  an  initial  slope  of  zero  and  then 
applies  a  linear  correction  to  the  deflections  to 
arrive  at  the  correct  deflection  at  each  nodal 
point . 

RUNG A :  This  subroutine  is  based  on  a  fourth  order  Runge- 
Kutta  method  of  integration  for  a  second  order 
differential  equation.  The  subroutine  calls  function 
"QF"  to  compute  shear  flow.  The  change  in  axial 
force  is  computed  between  every  pair  of  nodal  points 
by  integrating  the  shear  flow. 

QF :  This  function  defines  the  load-slip  relationship  for 

the  shear  connector. 

DEFLEX:  This  subroutine  computes  deflection  at  every  nodal 

point  by  using  a  second  order  Runge-Kutta  method  of 
integration  for  a  second  order  differential  equation. 

SYMSOL:  This  subroutine  is  an  equation  solver.  It  solves 

the  flexibility  matrix  to  compute  changes  in  reactions 
at  the  interior  supports. 

All  notations  for  main  vectors,  geometric  properties  and 

indices  are  defined  in  the  listing. 


START  ) 


READ  INPUT  DATA 


COMPUTE  SECTION  PROPERTIES 
(PROP) 


t 

COMPUTE  STRESS  RESULTANTS 
FOR  UNIT  REACTIONS 

AND  EXTERNAL  LOADING 
(CONC) 

* 

SUPERIMPOSE  STRESS  RESULTANT 
FOR  EXTERNAL  LOADS  AND 
SUPPORT  REACTIONS 
(CONC,  BMSU) 

f 

COMPUTE  SLIP,  SHEAR  FLOW, 
FORCE,  SLIP  STRAIN  AND 
DEFORMATION  FOR  EACH  UNIT 
REACTION  AND  STORE  THE 
FLEXIBILITY  COEFFICIENTS 
(MASTER) 

COMPUTE  SLIP,  SHEAR  FLOW, 
FORCE,  SLIP  STRAIN,  AND 
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